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Abstract 

We use the functional renormalization group equation for quantum gravity to construct a non- 
perturbative flow equation for modified gravity theories of the form S = J d d x^gf(R). Based on 
this equation we show that certain gravitational interactions monomials can be consistently de- 
coupled from the renormalization group (RG) flow and reproduce recent results on the asymptotic 
safety conjecture. The non-perturbative RG flow of non-local extensions of the Einstein-Hilbert 
truncation including fd d x^/gln(R) and Jd d x^/gR~ n interactions is investigated in detail. The 
inclusion of such interactions resolves the infrared singularities plaguing the RG trajectories with 
positive cosmological constant in previous truncations. In particular, in some i? _n -truncations 
all physical trajectories emanate from a Non-Gaussian (UV) fixed point and are well-defined on 
all RG scales. The RG flow of the ln(i?)-truncation contains an infrared attractor which drives a 
positive cosmological constant to zero dynamically. 
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1 Introduction 



Theories of /(i?) -gravity, in which the gravitational Lagrangian is based on an arbitrary func- 
tion of the curvature scalar, have recently attracted considerable attention within the cosmology 
literature. One of the prime reasons for this originates from the observation pQ that non-local 
terms which become dominant as R decreases can provide a natural explanation for the observed 
late-time acceleration of our universe, without the need of introducing dark matter (see, e.g., 
[2j for a review and more references). By combining these non-local interactions with higher 
derivative curvature terms, /(-R)-gravity can furthermore generate an inflationary phase in the 
early universe while at the same time satisfying the experimental bounds of gravity tests at solar 
system scales [3]. 

These phenomenological successes make it interesting to investigate /(i?)-gravity also from a 
quantum gravity perspective. A framework which naturally lends itself to such an analysis is the 
functional renormalization group equation (FRGE) for quantum gravity [4 J This equation allows 
a detailed investigation of the non-perturbative renormalization group (RG) flow of gravitational 
theories ^ |Zl 151 151 LH3 UU U21 US US1 US1 HZl H51 QSl I2D1 12S 122] , axid lias also been used to 
study the consequences of scale-dependent (running) gravitational coupling constants to black 
hole physics [23 [21], cosmology [23 [23 [23 [23 [22] , galaxy rotation curves [33 EE], and collider 
signatures of TeV-scale gravity models at the LHC [33 [33 131] • 

The main ingredient of the FRGE is the effective average action [35, 36J, a Wilson-type, 
coarse grained effective action which is constructed along similar lines to the ordinary effective 
action. Its construction starts from the standard generating functional for connected Greens 
functions W, obtained via the path integral. The usual action is then augmented by the addition 
of an infrared (IR) cutoff A^S", which provides a momentum dependent mass-squared term 7Zk(p 2 ) 
for the quantum field mode with momentum p: for p 2 3> k 2 , the cutoff function TZ^ip 2 ) vanishes, 
so that the high- momentum modes are integrated out in the usual way. For p 2 <C k 2 , we have 
lZk{p 2 ) oc k 2 so that the integration over the low-momentum modes is effectively suppressed by 
a mass-squared term proportional to k 2 . This modification induces a scale-dependence in the 
generating functional W^, which is then related to by a (modified) Legendre transform. 

One of the main virtues of Tf~ is that its scale-dependence is governed by an exact FRGE. 
Schematically, this equation takes the form 

ftr* = ±T* (1.1) 
2 r {2 > + n k 

with r^, denoting the second variation of Tj, and t := ln(/c). This equation defines a one- 
parameter family of effective field theories {T k ,k < k < k}, with k a UV cutoff scale and k, the 
effective or coarse graining scale. This family interpolates between the classical action S, imposed 
as an initial condition Tt = S at the scale k, and the coarse grained effective action which 
provides an effective description of physical processes with typical momentum k. Presupposing 



For a recent pedagogical introduction see [5]. 
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that the limit k — > exists (which is a priori not clear for a theory involving massless degrees of 
freedom), the standard effective action is obtained as T = Tq = lim^oT^. 

The framework outlined above provides a very powerful tool for quantizing a fundamental 
theory. In this case, one imposes the initial condition = S for the classical action S, uses 
the FRGE to compute for k < k and then takes the limits k — > and k — > oo. The defining 
property of a fundamental theory, loosely speaking, is that the limit k — > oo exists. Typically, 
taking this limit requires the existence of a UV fixed point of the RG flow. In the context of the 
FRGE for quantum gravity, there is by now good evidence that such a Non-Gaussian fixed point 
(NGFP) indeed exists, defining a predictive fundamental theory of gravity along the lines of the 
asymptotic safety scenario proposed by Weinberg [37] (see [381 ESI SQj for recent reviews) §|. 

Moreover, the FRGE can also be used to investigate the structure of the gravitational effective 
action in the IR irrespectively of the existence of a fundamental theory. Here, one starts from 
an effective theory at the scale k and uses the FRGE to evolve the theory down to scales k < k, 
successively integrating out the quantum fluctuations "shell by shell" . In particular, it allows one 
to address the proposals made in [45\ l4Uj |4"7] that quantum gravity gives rise to strong RG effects 
in the IR which are responsible for the tiny value of A. At the level of the gravitational RG flow, 
this would most likely be reflected by the existence of an IR fixed point, which attracts the RG 
flow for k — * 0. Such an IR fixed point could provide a dynamical solution to the cosmological 
constant problem, if A is driven to zero irrespectively of its initial value at the scale k. 

The main limitation of the FRGE, however, is that it cannot be solved exactly. In order to 
extract some non-perturbative information from it, one therefore resorts to the approximation 
technique of truncating T^. That is, one makes an ansatz for which only retains a finite (or 
infinite) subset of all possible interaction monomials compatible with the symmetries of the theory. 
Introducing ^-dependent coupling constants multiplying the interaction monomials, this ansatz is 
substituted into the exact FRGE. By projecting the resulting equation onto the subspace spanned 
by the truncation, one can then construct the /3-functions for the scale-dependent couplings 
without evoking perturbation theory in a small parameter. Depending on the number of couplings 
included in the ansatz, this procedure leads either to a coupled system of ordinary differential 
equations (finite number of couplings) or to a partial differential equation (infinite number of 
couplings) governing the scale-dependence of the truncated r&. Consistency of the truncation 
ansatz thereby requires the extraction of all terms proportional to the interaction monomials 
spanning the truncation ansatz from the traces appearing on the RHS of the FRGE. 

Finding good truncation subspaces generally poses a delicate problem, requiring some physics 
intuition about the interaction terms which are relevant for describing the physics at momentum 
scale k 2 , and usually involves some educated guesswork. Contrary to perturbation theory, where 
the quality of the approximation is controlled by an expansion parameter, approximating the full 

2 Some features associated with the NGFP have also been observed in the apparently orthogonal Causal Dy- 
namical Triangulations (CDT) approach to quantum gravity [41 II42||4"5] . thus possibly hinting at a hidden relation 
between the FRGE and CDT. Moreover, the existence of the NGFP is also supported by studies of the symmetry 
reduced Euclidean path integral [44] . 
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Figure 1: The RG flow of 4-dimensional quantum gravity in the Einstein-Hilbert truncation. 
The arrows point in the direction of increasing coarse graining, i.e., towards decreasing k. 

RG flow of the theory by a truncation does not allow one to judge the quality of the truncation a 
priori§| So far, most work on gravitational RG flows has been carried out in the Einstein-Hilbert 
truncation, in which the gravitational part of Tk[g] is of the Einstein-Hilbert form, 



with Gk and denoting the running Newton's constant and cosmological constant, respectively. 
The associated /3-functions governing the RG flow of the corresponding dimensionless couplings 
9k = Gkk d ~ 2 , Afc = Afc//c 2 , have been first derived in [3], and its solutions, classified in [8], 
are shown in Figure [TJ There, the most prominent features are the UV attractive NGFP at 
g* > 0, A* > and the boundary of the coupling constant space at A = 1/2, which causes the 
termination of the type Ilia trajectories (positive cosmological constant in the IR) at a finite value 
&tcrm > 0. Now, due to the truncation approximation, there is no a priori reason that the NGFP 
and the A = 1/2-line carry over into the full theory space. Indeed, while there exists substantial 
evidence that the NGFP is not a truncation artefact and is present in the full theory, we will 
show that the A = 1/2- line can be generically resolved with the inclusion of non-local curvature 
terms in the truncation ansatz. 

One of the main results of this paper is the extension of the /3-functions obtained in the 
Einstein-Hilbert truncation to a partial differential equation governing the non-perturbative RG 

3 A posteriori, however, there are indirect tests, such as, e.g., analyzing the (unphysical) cutoff function depen- 
dence of observable quantities [5]- 




(1.2) 
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flow of /(-R)-gravity theories of the form 

T k [g] = (IQnGk)- 1 j d d x^f k (R) , (1.3) 

where fk(R) is an arbitrary, scale-dependent function of the scalar curvatureQ This is achieved by 
performing the transverse-traceless (TT) decomposition of the metric fluctuations and combining 
the construction of the FRGE with "type B cutoff" [7] (which we will denote as TT-cutoff in the 
sequel) with the ideas outlined in [22\ . Subsequently, we use this flow equation to study the RG 
flow of specific non-local extensions of the Einstein-Hilbert action 

f k [g] = J^r J d d x^fg-{-R + 2A fc + 16vrG fc v k ln(i?)} , (1.4) 

and 

f k[g] = Y^^jd d x^{-R + 2A k + 16nG k fl k R- n }, n>l€Z. (1.5) 

The non-local terms included in the truncations (|1.4j) and 111. 5ft serve as prototypes for interactions 
which become dominant for small values of the curvature scalar. In this respect, the primary goal 
of this paper is the qualitative understanding of the RG flow for the new non-local couplings v k ,fi k - 
The truncations (II. 4p and (jl.5j) thereby constitute simple models which allow one to investigate if 
and how non-local interactions of this type affect the gravitational RG flow. Demanding that the 
truncated T k give rise to stable cosmological solutions, which (most likely) do not exist for classical 
actions of the form ()1.4[) and (jl.5p . requires the extension of the truncation ansatz by including 
higher (positive) powers of the scalar curvature. Based on the work on i? 2 -truncations we expect 
that these higher derivative couplings can appear rather naturally along the RG trajectory, so 
that passing on to phenomenological acceptable models is straightforward once the dynamical 
origin of the non-local interactions is established. 

Studing the RG flow based on the truncations (|1.4p and (|1.5p also constitutes a natural 
extension of the results |13j . where theories containing non-local volume terms were considered. 
The study of the RG flow of the couplings v k ,ji k hence provides a natural next step in this 
program, laying the groundwork for analyzing truncations including the (probably more physical) 
non-local interactions f d d x^/gRln(—D 2 )R and J d d x^/g~R(—D 2 )~ l R proposed in [45] and [50] . 

The rest of the paper is organized as follows. Section [2] constructs the FRGE for quantum 
gravity utilizing an improved geometric gauge-fixing scheme and transverse traceless decompo- 
sition. This FRGE is restricted to /(i?)-gravity in Section and the tools for evaluating the 
corresponding operator traces are there introduced. Section [4] discusses the possibilities for de- 
coupling certain gravitational interaction monomials from the RG flow and the resolution of the 
IR singularities appearing in the Einstein-Hilbert truncation. Section [5] and [6] contain the detailed 

4 At the classical level (or for fixed values k), certain gravitational actions of this form can be mapped to the 
Einstein Hilbert action coupled to one additional scalar field [48] . While the later formulation is very useful in 
analyzing the solutions of the theory, the conformal rescaling which relates the two dual formulations introduces 
a fc-dependence of the metric and scalar field in the dual frame, making the derivation of the RG equation in the 
dual frame rather cumbersome. Thus, we will derive the RG equation based on fk in the following. 
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analysis of the RG flow arising in the truncations (11.41) and (jl.5p . respectively. Lastly, using our 
results from Section [31 we evaluate the operator traces for a special choice of cutoff function in 
d = 4, deriving a partial differential equation governing the RG flow of /(i?)-gravity in Section 
before concluding with a short summary and discussion of our results in Section [8l The technical 
details on the heat-kernel techniques and threshold functions employed in the construction of the 
RG equations are collected in Appendix [A] and [Bl respectively. 

2 The improved FRGE for quantum gravity 

We start by constructing the flow equation for /(i?)-gravity, rederiving the FRGE for quantum 
gravity based on the TT-cutoff [7] and including the improvements [22]. The novel features of 
our construction are the implementation of a distinguished, geometrically motivated gauge-fixing 
condition and the sidestepping of the momentum dependent field redefinitions appearing in the 
original construction. As pointed out in [22], this leads to a significant simplification of the 
TT-cutoff. 

2.1 The geometric gauge-fixing condition 

In order to assure that our FRGE respects diffeomorphism invariance, we follow [3] and define 
a scale-dependent Euclidean functional integral using the background field method. Thereby, 
the quantum metric 7^ is decomposed into an arbitrary but fixed background metric and a 
(not necessarily small) fluctuation field h^ u 

7fj,u = g^u + . (2.1) 

The formal expression for the generating functional then includes integrations over the quan- 
tum fluctuations h^ u and the Faddeev- Popov ghosts C^, C^, 

Z fc [sources] = / Vh^VC^VCn 

J (2.2) 

x exp [S[h + g] - S gi [h;g] - S gh [h,C,C;g] - A k S[h,C,C;g] - S SOUTCe ] . 

Here, S[j] = S[h + g] is a diffeomorphism invariant action functional which, for the time being, 
is assumed to be positive definite. Furthermore, S g f[h;g] denotes the gauge-fixing term, 

% l h ->9] = ^J d d x^ r v ^ [h; g] F v [h; g] , (2.3) 

which depends on the gauge-fixing parameter a and implements the gauge-fixing condition 
F^h, g] = 0. For the construction of the FRGE, it is particularly convenient to choose a F^[h; g] 
which is linear in h^ u . Explicitly, we take 

F m [/i;#]=V2k(^V-^A^V)> ( 2 - 4 ) 

with k = (327rG) _1 / 2 . Here and in the following, barred and unbarred quantities are constructed 
with respect to the background metric and the full quantum metric, respectively. In particular, 
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D p denotes the covariant derivative based on g pu , while D p is covariant with respect to "f pu . The 
constant p parameterizes the freedom of implementing different gauge choices, with the harmonic 
or de Donder gauge employed in [7j being obtained by setting p = | — 1. In the following, it will 
however turn out that it is more convenient to set p = 0, as this results in drastic simplifications. 
The ghost action arising from the gauge fixing (|2.4p is constructed in the standard way and reads 

S gh [h,C,C;g] = -V2 j d d x^C ti M" u C u , (2.5) 

with the Faddeev-Popov operator 

M\ = Dp rf v D p + DP lpv DP - 2 D p gP° lpv D a . (2.6) 

The remaining terms in (|2.2|) . Af.S[h, C, C; g] and 5g OU rce) contain the scale-dependent IR cutoff 
and the source terms for the metric fluctuations and ghosts, respectively, and will be discussed 
in Subsection 12.31 

2.2 The transverse traceless decomposition 

One of the main technical difficulties in practical applications of the FRGE for gauge theories and 
gravity is the construction of the IR cutoff AkS[h, C,C;g] and the inversion of rjf + K k , which 
requires the (partial) diagonalization of these operators. With the TT-cutoff, this is facilitated 
by performing the transverse-traceless (TT) decomposition [51] of the metric fluctuations and the 
ghost fields^ Here, the metric fluctuations h pv are decomposed into a trace part h 1 * = ^g^cf), 
encoded in the scalar cj>, a transverse vector the longitudinal part of which is given by a, and 
the transverse traceless tensor K^ v according to 

hpv = h pu + D^v + D u ^ + D p D v a - ^g pv D 2 a + , (2.7) 

with the component fields subject to the constraints 

g^h^ = 0, S"/£, = 0, S%, = 0, tf> = g (lv h^. (2.8) 

Similarly, the ghost fields C^,C M are decomposed into their transverse and longitudinal parts 
C T ",Cj and rj, 77 as 

C fi = C$ + Dtf, C p = C^ + D pV , DPCJ = 0, DPCj = 0. (2.9) 

From eqs. (|2.7p and (|2.9p . it is obvious that not all modes of the component fields contribute 
to the metric and ghost fields. In the metric decomposition, the constant mode of the cr-vectors, 
C p = D p a satisfying the conformal Killing equation 

D p c u + b v Cp. - - d 9^b^ = , (2.10) 

5 This decomposition is optional when considering the Einstein-Hilbert truncation with a = 1, but mandatory 
when higher powers of the curvature scalar are included in the truncation ansatz. 
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and transversal vectors solving the Killing equation 







(2.11) 



do not contribute to h^ u . Analogously, the ghost fields C^,C^ are independent of the constant 
modes of rj, fj. These modes are unphysical and must be excluded from the spectrum. In fact, this 
is also a necessary requirement for the invertibility of the Jacobians arising from the coordinate 
transformations (|2.T[) and (|2,9p . as we will see below. 

A virtue of the TT-decomposition is that the component fields (almost) provide an orthogonal 
basis for the quantum fluctuations. Evaluating the scalar product in the ghost sector gives 

(C,C) := f d d x^gC^ = f d d x^{ClC^ -fjD 2 V }, (2.12) 

establishing that this decomposition is orthogonal. In the metric sector, on the other hand, the 
inner product becomes 

( h W,hW):= fd*xy/§h$g"p>h@ 



d d x^[h$ T h^ T ^ - 2£« [g^D 2 + iff") £( 2 ) - 2£«i™^ {2) (2-13) 

- 2^R^D u a^ + a® [^-(D 2 ) 2 + D^D u ]a^ + ^ (1 ¥ 2) ] , 

which is orthogonal up to the ^-cr-mixed terms. We note that these crossterms vanish, how- 
ever, for the Einstein backgrounds (satisfying R^ u = Cg^ u ) considered in the next section, 
/ d d Xy/g^W v D v a = -C f d d x y /gaD' M ^ = 0. Thus, in that case, the component fields do 
provide an orthogonal basis for the quantum fluctuations. 

The inner products (|2.13p and (|2,12p can be used to compute the Jacobians of the coordinate 
transformations (|2,7p and (|2,9p . For this purpose, we follow |52] and consider the Gaussian 
integral over and C^, C M . In the metric sector, this yields 



Vhp„ exp [~l(h, h)] 



J grav y VhlPZyPoVif, exp [ - \ j d d x^g {fcJ> T "" + ±0 2 + a] [&,, a] T ] 



(2.14) 



with M^' u ^ a (d+ 1) x (d+ l)-matrix differential operator whose first d columns act on transverse 
spin one fields £u and whose last column acts on the spin zero fields a. The corresponding matrix 
can be read off from (|2.13p . 



-2 [g^D 2 + BP"] 



-2R>* X D> 



2D\ R Xv ^{D 2 ) 2 + D^R^D, 



(2.15) 



The Jacobian J grscv is found by performing the Gaussian integrals in (|2.14p and, up to an infinite 
normalization constant, which can be absorbed into the normalization of the measure, reads 



J, 



grav 



det i 



(1T,0) 



1/2 



(2.16) 
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Here, the subscript (IT, 0) refers to the matrix structure explained above and the prime indicates 
that the unphysical modes are left out. The Jacobian introduced by the change of coordinates in 
the ghost sector ()2.9p can be determined analogously. Setting 



PC^PC M exp[-(C,C)] 



Jgh JvC T ^VC^Vr,Vfjexp - J d d x^g {Cj - fj D 2 v } 



(2.17) 



and performing the Grassmann functional integrals results in 

J gh =(det' [-D 2 \y 1 . 



(2.18) 



In order to include the contribution from the Jacobians in the partition function (|2.2p . we 
essentially copy the Fadeev-Popov trick and introduce a set of auxiliary fields to exponentiate 
the determinants. This leads to an additional term Sauxfaux. fields; <?] in Let us first consider 
Jgrav arising from the metric sector. Introducing the transverse ghost cJ^,c TfJ- , a "longitudinal" 
Grassmann scalar b, b, a transverse vector field £j and a real auxiliary scalar to, we can use 



standard results on Gaussian integration to write (|2. 16|) as 

Jgrav = J V^Vc^^VbVbVQVu 



x exp 



d d xV=9 {[C>] [MM]'[gM T + Si [M M ]'[cl bf} 



(2.19) 



The Jacobian from the ghost sector (|2.18p is exponentiated by introducing a new complex scalar 
field s, s, 



Jgh = J exp 
By combining the results (|2.19p and (|2.20p . we then find 

5 , aux [C T ) s, s,uj,c t ,c T ,b, b;g] = 



d d x^s[-D 2 }' s 



(2.20) 



(2.21) 



Here, the primes are placed as a reminder that the unphysical modes have been excluded from 
the auxiliary fields. 



2.3 The FRGE in terms of component fields 

After discussing the new ingredients for the improved FRGE with TT-cutoff, we are now in 
the position to construct the corresponding flow equation. To lighten our notation, let us first 
introduce the index sets 

h = {h T ,Z, a A} , h = {C T ,C,fj, V } , I 3 = {( T ,u J ,s,s,c T ,c T ,b,b} 1 (2.22) 
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and a shorthand notation for the quantum fields, 

X = {/i T ,e^>,C T ,^s, S ,C,77, C T ,6,C T ,f?,c T j}. (2.23) 



The two remaining ingredients in the scale-dependent partition function (|2.2p which have not 
yet been specified are the IR cutoff A^S and the source terms. The IR cutoff is constructed in such 
a way that, for all fields (including the auxiliaries), the integration over the p 2 = — l) 2 -eigenmodes 
with large eigenvalues p 2 3> k 2 are unaffected while the contributions of the modes with small 
eigenvalues p 2 <C k 2 are suppressed. To this end, A k S provides a momentum-dependent mass 
term which is quadratic in the fields (|2.23p 

A k S = ±Yl (Ci,(^f av )c 1 c 2 C2) + i £ <V>i,(^f h WV>2> + ± £ (a.WTWa)- ( 2 - 24 ) 

CliC2e/l 1pl,t(>2£h <fl,«Gi"3 

The operators 7^f rav , and 7?.^ ux depend on the background metric only. In order to implement 
the desired suppression behavior, they must vanish for p 2 /k 2 — > oo (and in particular for fc 2 — > 0) 
while, for p 2 /k 2 — > 0, they must behave as !Z k - > Z k k 2 , with Z k a possibly matrix valued In- 
dependent function. In addition, hermiticity requires that (7Zk)( 2 (i = C^fc)^^ an< ^ C^-k)ip 2 i>i = 
— (^fc)L^ 2 f° r bosonic and Grassmann valued fields, respectively. Furthermore, (7Zk)ipnp 2 = 
if both ip\,ip2 G {C T > 7 7> cT >fr} or {C lT > 77, c T , 6}. In the sequel we will denote the three terms in 
(|2.24p by AfcS'grav, AfcS'gh and A fc S aux , respectively. 

Finally, we specify Ssource; which introduces source terms for all component fields, again 
including the auxiliaries. Denoting the sources for the bosonic field and ghost fields by J and K, 
we set 



Source = ~ ^(J^X) ~ (^C T '£ T ) ~ (Ju),^) ~ (Js,s) - (J s ,s) 

<eii 

^e{C T ,»7,c T ,fe} VG{C" r ,r?,c T ,f)} 



(2.25) 



We hereafter denote the sources collectively by 

J = \^Jh^^^^o^<i>iJc^ u '^^ s ^c T ^v^c T ^ h ,K^T,K n ,K Z T 1 K^ , (2.26) 
which allows us to write (|2.25|) schematically as 

S SO uvce[J;9] = -(J,X)- (2-27) 

The relations between the component sources and the standard sources for the metric fluctuations 
h^u, C M have been worked out in detail in [7], where it was shown that the standard sources 
can be reconstructed from (|2.26p . 

With these results at hand, we can now write the scale-dependent partition function (|2.2p in 
terms of the component fields, 

Z k [J;g} = I Vh T V^VaV(f)VC T VujVsVsVC Tll VC^VriVfjVc T ^VclVbVb 

J (2.28) 

x exp [-S[h + g] - S gi [h; g] - S gh - A k S - 5 aux - S sourcc ] . 
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Note that the fc-dependence on the RHS of this equation arises solely through the IR cutoff A k S. 
The derivation of the improved FRGE then proceeds completely analogously to [7\. Based on 
(|2.28p . the scale-dependent generating functional for the connected Green's functions is given by 

W k [J; g]= In Z k [J;g\. (2.29) 

From this, the classical component fields are obtained as 

/ \ 1 5Wk rt nM 

W - (Xi) = Tg W , (2-30) 
where we denote the classical counterparts of the quantum fields (|2.23p collectively by 

{^ T , i, v, 4>, C T , w, t, t, v T , g, u T , t, v T , g, u T , f | . (2.31) 



The classical counterparts of the fundamental fields h^, and can be reconstructed by 
substituting the corresponding component fields into the formulas for the TT-decomposition. 
For later purposes, we also define the classical analogue of the quantum metric (|2.ip 



(2.32) 



Subsequently, the scale-dependent effective action is constructed as the Legendre-transform on 
W k with respect to the component sources J, taking the latter as functions of the classical fieldfi 



F k [<p;g] = (J,<p)-W k [J;g] 



(2.33) 



The effective average action T k is then defined as the difference between and the cutoff action 
with the classical fields inserted, 



Fkiwg] '■= ffc^ff] - AkS[<p;g] . 



(2.34) 



The derivation of the FRGE governing the fc-dependence of r^f^g] proceeds in several steps. 
First, we differentiate (|2.29|) with respect to t = ln(k) 



-d t W k = ±Tr' 



{ipi ® -02)<9t(7?-fc)v>ii/>2 



+ ±Tr' 



Y (Xi ®X2)dt{K k ) Xx 

Xi,X2&h 



x> 



(2.35) 



where we expressed the RHS in a matrix notation. We then introduce the Hessian of the effective 
action, 



r£ 2) )>,y) := (-i) 



1 



(2.36) 



6 As was shown in [7J, taking the Legendre transform with respect to the sources for the (physical) component 
fields is equivalent to taking the Legendre transform with respect to the sources for the fundamental fields. 
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with [j] = 0, 1 for commuting fields (fj and Grassmann fields (fj, respectively. This Hessian is the 
inverse of the connected two-point function 

(G k)ij ( X ,y) := ( XiXj ) - nWvM = ^^ SjWm ' (2 ' 37) 
in the sense that 

'dfiyy/M (Gkhj (x,y) (f^Y (y,z) = 5<^=^ . (2.38) 

Using these relations, we can express the expectation values (Xi( x )Xj(y)) i n ec l- <|2.35j) through 
(f^ ) (y> z ) + 4 > i{ x ) ( t ) j{y)- Finally, performing the Legendre transform f)2 .33j) and subtracting 
AkS[(p;g] yields the desired FRGE for r^f^g] 



+ ^ T 4 E (r^+nX 1 d t (n k ) M . 



(2.39) 



L - V ' X1X2 

Here, we again used a matrix notation on the RHS, and the barred index sets run over the classical 
fields, 

h = {/l T ,e>,0} , I 2 = {V T ,V T ,Q, Q } , I 3 = {( T ,LU,i,t,U T ,U T ,f,T}. (2.40) 

The first, second and third term in (|2.39[) thus encode the contributions from the gravitational, 
ghost and auxiliary fields, respectively. 

Eq. (|2.39p is the desired exact FRGE for the effective average action employing the improved 
TT-cutoff. Compared to [7], the new feature is the appearance of the third trace term, capturing 
the contribution of the auxiliary fields, which replaces the momentum dependent field redefinitions 
in the metric and ghost sector. As it will turn out in the next section, this form of the FRGE is 
very convenient when studying truncations involving a general function f k (R). 



3 Constructing the flow equation of /(i^)-gravity 

After deriving the exact FRGE (|2,39p in the last section, we now proceed by projecting the 
resulting RG flow on truncation subspaces spanned by arbitrary functions of the curvature scalar. 
In view of recent applications of RG methods to gravity theories including extra dimensions, we 
do not fix the space-time dimension d and keep the shape of the IR cutoff generic. 

3.1 The truncation ansatz 

Let us start by specifying our truncation ansatz for T k [ip; g], which we take to be of the form 

T k [ip;g] = T k [g] + S gi [g -g;g} + S gh [g - g,C T ,C T ,fj,rj;g] + 5 aux , (3.1) 
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with 

r k [g] = Y^G~J ddx V9fk(R) ■ (3.2) 

This ansatz captures the RG flow of /(i?)-gravity by promoting the classical function f(R) to be 
RG scale-dependent. The RG effects in the gauge-fixing, ghost and auxiliary sector are neglected 
by treating 5gf,5 g h and 5 aux as classical, i.e., not containing any scale-dependent couplings. A 
similar form of this truncation ansatz (without resorting to the auxiliary field formulation) has 
already been considered in the pioneering paper [3J, where it was shown that neglecting the RG 
running in Sgf and S g h is, to leading order, compatible with the modified Ward identities of the 
theory. Moreover, even when treating SW as classical, the RG effects in this sector can be included 
by taking a = 0, which corresponds to a fixed point of the RG flow |53| . 

When substituting the truncation ansatz (|3.ip into the FRGE (|2.39[) we find that the gravi- 
tational trace only receives contributions from 

rr v [g;g}:=T k [g]+S si [g-g;g], (3.3) 

while the ghost and auxiliary traces are evaluated based on the classical actions S g h and S anx , 

d t T k [^g]=\Tr>[ £ (rf w(2) + K k ) ~\ d t (K k ) ClC 



C1.C26/1 



+ §M E (^t + KkY 1 d t {K k ) XlX2 



(3.4) 



X1.X2G/3 



The index sets here are given by (|2,40p . 

Our main task now is to evaluate the traces appearing on the RHS of (|3,4p . This is done as 
follows. We first compute the Hessians of the various terms appearing in (|3.ip in Subsection 13.21 
This computation can be simplified by choosing the background metric g^ as the one-parameter 
family of metrics on the (i-sphere S d (parameterized by the radius r or, equivalently, the curvature 
scalar of the sphere), as this suffices to distinguish different functions of the curvature scalar. For 
these backgrounds, the unphysical modes discussed in Subsection 12.21 correspond to the lowest 
— Z) 2 -eigenmode (I = 1) of the transverse vector the two lowest eigenmodes (I = 0, 1) for 
a, and, similarly, the lowest mode (I = 1) of the transverse ghost field C T/i , Cj as well as the 
two lowest modes (I = 0, 1) for the scalar ghost 77, f]. Using Table H] in Appendix [A] one can 
explicitly check that these modes are annihilated by the operators appearing in the Jacobians 
(I2TTU1) and (ETTSfl restricted to S d . In Subsection E3J these results are used to construct the IR 
cutoff operators AS k explicitly. Substituting those results back into (|3.4p . we find cancellations 
between several terms appearing on the RHS, so that the final result (|3.2ip takes a surprisingly 
simple form. 
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(2) 

3.2 Computing the Hessian 1^ 

We start by deriving the term quadratic in the metric fluctuations h arising from Tk[g = g + h;g], 
(|3.2p . For this purpose, we expand T k [g + h;g] in a Taylor series in h, 



T k [g + h;g} = T k [g; g] + 0(h) + T^ A [g + h-g]+0(h i ). 



(3.5) 



Here, T^ uad [g + h;g] is found by taking the second variation of T k [g + h;g] with respect to h and 
setting g = g afterwards. Applying the chain rule, we obtain 

1 



s 2 t h 



167rG fc 



d d x 



f k (R)5 2 ^ + f' k (R) {2(6Jg)(6R) + y/g{PR)) + ^gf' k \R){5R) 2 , (3.6) 



where, here and in the following, the prime denotes a derivative with respect to R, i.e., f' k (R) = 
9 ^8R ' e ^ c ' ^ e then substitute in the above expression the variations of the metric and Ricci 



scalar, which can, e.g., be found in [54] . Note that the result drastically simplifies once we set 
g = g with a spherically symmetric background metric. Performing the TT decomposition of L 
according to (|2.7p . a lengthy but straightforward computation yields 



=,quad 



[h;g] 



with 



2 fW 



Jfe ^ - Jk + S(d^T) u Jk 



+ i (df k -2Rf' k ) ^ 



D 2 + \R 



32nG k 



^- J d d x^C k , 



(3.7) 



+ 



Ad- 



4(d - l) 2 /^ 4 - 2(d - 1) ((d - 2)# - 4 5 /£) D 2 



(3.8) 



+ (d-2)(d/ fc -4i^)+4 J R 2 //' 
2(d - lff^D 8 - (d - l)((d - 2)/£ - 4i?/£)L> 6 

- (d(d - 1)/* - + 2Rf£)) D 4 - R(df k - 2Rf' k )D 2 

1 - 2{d - \ffl& -{d- l)((d - 2) # - 4tf/£)L> 4 - fl((d - 2)f k - 2Rf k >)D 2 



a 



"d 1 



a . 



Here the barred quantities are constructed from the background metric g and we have suppressed 
the argument of f k (R). 

Following the procedure outlined above, we also extract the quadratic terms arising from 
<Sgf , 5* g h and Saux- Applying the TT decomposition in the first two cases yields 



quad 
gf 



K 

a 



d d x 



D 2 + -R 



d 



((d- i)d 2 + r) 2 jy 



(7 



(3.9) 



P_ 

"d 2 



quad 
gh 



(d- l)D A + RD 2 a}, 
-V2 J d d x Vd{ [D 2 + l d R] C ^-^v[(d-l-p) D 4 + RD 2 ] V ] , (3.10) 
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while substituting the spherical symmetric background metric simplifies S anx to 

^x ad = / d*Zyfi{C T *[-&-$R]\* + 

J L (3.11) 

+ s [-D 2 }> s + b[D A + ^RD 2 ] "b + u;[D A + ^RD 2 ] " u) . 

Here, the primes on the operators indicate that the corresponding number of lowest eigenmodes 
are excluded from the physical fields. Note that, for the special case d = 4, eqs. (|3.7p . (|3.9p 
and (|3.10p are in precise agreement with earlier findings [55J. We also remark that, while in 
the harmonic gauge used in [7] performing the TT-decomposition in the ghost sector is optional, 
working in the gauge (|2.4I) requires this decomposition in order to diagonalize the quadratic 
fluctuations. 



3.3 Adapting the cutoff operators 

Our next task is to explicitly construct the IR cutoff operators ()2.24p . For this, it is useful to 
split 1Zk{p 2 ) into its matrix part and a scalar function Rj~, encoding the momentum-dependent 
mass term, 

[**(PX W = i Z *U„ R ^ = t^U 2 kVH-D 2 /k 2 ) . (3.12) 
The dimensionless profile functions R^°\—D 2 /k 2 ) interpolate between R(°\Q) = 1 and R(°\oo) = 
and otherwise allow one to adjust the "shape" of the momentum-dependent mass-squared term 
- see App endix IB . 2 1 for an explicit choice of a profile function specifying a cutoff-scheme. The aim 
of the IR cutoff is the regularization of Hessians 

[^k\ lV2 = [f(-D 2 ^...)] viV2 (3.13) 

in such a way that all covariant derivatives appear in the form 

P k = -D 2 + R k (p 2 ) . (3.14) 

Thus, as a result of the IR regularization, 

l?f +M^ = [ft- D2 + (3.15) 

depends on the covariant Laplacian — D 2 through the combination (|3.14p only. 

With these prerequisites, we can now explicitly construct A k S for the truncation (|3.ip . In 
order to simplify the resulting expressions, we use the gauge- freedom retained in (|2.3p to set0 

p = 0, a = 0. (3.16) 

As eq. (|3.19p below will illustrate, this gauge choice leads to a factorization of the physical and 
gauge degrees of freedom contained in h^ u in such a way that the physical degrees of freedom 

7 We note that the limit a — > has to be implemented with care and, strictly speaking, should be taken only 
after constructing the flow equation to avoid singular expressions. In any case, a more careful treatment of this 
matter confirms the validity of the final result (|3.2ip . 
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are parameterized by h^ u and <fi, while the gauge degrees of freedom are confined to £j and 
its longitudinal component, a. Thus, from a geometrical point of view, f)3. 16j) aligns the gauge 
orbits with the directions £j,cr. This feature will later lead to a considerable simplification of 

(2) 

the flow equation. In particular, we will find that the matrix structure of the operator T, + lZ k 
diagonalizes, so that the latter can be easily inverted. 

This said, let us start with constructing the IR cutoff for the auxiliary fields. Considering 
(|3.11|) . it is straightforward to establish that the prescription (|3.15|) is implemented by 



A fe S aux = J d d x^g {c T/i i?feCj + c T »R k cl + sR k s + b [Pi - ^RRk - D 4 

■}■ 



+ C0 



(3.17) 



Analogously, starting from (|3.10j) . we find the IR cutoff in the ghost sector, 



A fc ,S gh = d d x^§ [d^RkC^ + 



2(d-l) - 



P 2 

I I. 



d-l 



RRi 



(3.18) 



The construction of A^/Sgrav is slightly more involved. Here, we first observe that, for p = 0, eq. 
(13.91) diagonalizes and contains the fields £ M and a only. Combining (13. 7ft and (13.9|) then gives 

1 



-igrav;quad 



327rGfc 



(TXy 



^{\hl v [f k &-f k + ^Rf k ]h^ 



(3.19) 



+ 



id- 



4(d - l) 2 /^ 4 - 2(d - 1) {(d - 2)f' k - 4 R /£) & 



+ ^ J d d x^g{^ 



+(d-2)(df k -4Rf k )+4R 2 f' ! !] cj>) 



D 2 + -R 

a 



((d-l)D 2 + R) 2 D 2 



a + O 



(«)} 



^quad 



d? 

Compared to the terms appearing in Sg, the ££,cr<7 and c/xr-contributions from are sup- 

pressed by one power of a and are indicated by 0(a) in the formula above. In the limit a — > 0, 
these terms do not contribute to the flow equation and will therefore be neglected in the following 
discussion. Again following the rule (|3,15p . we find the IR cutoff in the metric sector, 



1 



32vrGt 



d d x^\±h T 



2 V 



(3.20) 



+^ 4> [4(d - l) 2 f' k '{Pl - £ 4 ) + 2(d - 1) ((d - 2)f' k - ARf") R k ] cp} 



~ J d d x^r 9 {^[pl-D^-^RR k 



+ d^ a 



[(d - l)P k -R) 2 P k - {{d - 1){-D) 2 - R) 2 (-D 2 



a + O 



(«)}• 



This result completes the adaptation of the IR cutoff to the truncation (|3.1j) . From the results 
(|3.17|) . (|3.18|) and (|3.2U|) . it is then straightforward to read off the entries of the cutoff operator 

(2) 

lZ k and construct TV +7Z. k , which indeed assume a diagonal form in field space. For convenience, 
their entries are summarized in Table [TJ 
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Cutoff operator IZk 



Matrix element of rj? + IZk 



32-KG k f'k R k 
64TrG k d 2 /v fc 



167rGoa 



p2 



D 4 



y,RRk 



V2R k 

2R k 

Rk 

Rk 



[Ph 



p2 



^Rfp.+ ^ + ^RfD* 



D 4 



d=l RR k 



p2 



d-l 



RRi 



P k ~D A - d=T RR k 



l 

327rG k 



fkPk-fk + ^ R f' k 



i r 

mnG k d' 2 k 



(2)<t>4> 



167rGoa 

(d-l) 2 
16-KGoad 2 



hR 



{ R k - dzjRfPk 



V2 
2 



Pk~\ R 



Pi 



-dh R 



Pk-jiR 



h-hR 



Pk 
P k 



d^T R 

- dh R 



Pi- 



(2) 

Table 1: Matrix entries of the operators IZk and T k + 1Zk to leading order in a. The first column 
indicates the indices of the matrix element in field space, while the second and third columns 

(2) 

contain the corresponding matrix element of IZk and T k + IZk, respectively. The elements are 
symmetric under the change of bosonic indices, while they acquire a minus sign when Grassmann- 
valued indices are swapped, itff' and rj^ 9 ^ are defined in eqs. (|3.22[) and (|3.23[) . 



With these results, we can now evaluate the traces in (j3.4f) for the truncation (|3.ip . Substi- 

1 v 1-1 listed in Tableland properly taking into account the 



tuting the operators IZk and (r^ + IZk 
non-physical modes, we fin 



u k 



2 lr 



d t Rk 



Pi 



-j-^rR 



d-l 



1 



+ £>i(d,0) -^dtRk 

Pk 



-D 2 =Ai(d,0) 



l Tr / 
2 ±r lT 



d t Rk 



Pk 



\ R 



+ ^Tr 2T 



d t {Z Nk f k Rk) 



%Nk [flPk + fk 



2(d-2) 
d(d-l) 



R f'k 



+ iTr 



d t {z Nk izt) 



ZNkF 



{2)4><t> 
k 



(3.21) 



with f k and fk(R) introduced in eq. (13.21) . Here, Z^k captures the scale-dependence of Newton's 
constant via Gk ■= Zj^Gy with G k denoting its "bare" value. The subscripts 0,1T,2T on the 
operator-traces indicate that their arguments act on fields with spin s = 0, 1, 2. Furthermore, the 
primes imply that the corresponding number of lowest — Z) 2 -eigenvalues are excluded from the 
trace. Following the notation of Appendix |Aj the multiplicity and — Z) 2 -eigenvalues of a given 



8 A similar equation also underlies the analysis of the UV properties of /(-R)-gravity [22] . The corresponding 
equation differs from (|3.2HI by the contribution of the single eigenmode proportional to Di(d, 0). We thank R. 
Percacci and C. Rahmedc for discussion on this point. 
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field mode are denoted by Di(d, s) and Ai(d, s). For the d-sphere, these are summarized in Table 
[4] in Appendix [A] Finally, we have introduced 



TZf = 4 (d - lff k > (P 2 - (-D 2 ) 2 ) + 2 (d - 1) {(d - 2) f k - AR ff) R k 



(3.22) 



and 



ff** = 4 (d - lff k ' (p k - ^R) + 2 (d - l)(d - 2) f' k (p k - ^R) + d(d - 2)f k , (3.23) 



in order to write the last scalar trace in compact form. After carrying out the variations 
T k [g, . . . ;g], we can also identify g = g in (|3.2ip . dropping the bar on the metric in the sequel. 
For later reference, we denote the five terms appearing on the RHS of (|3.2ip by S l k , i = 1, . . . , 5, 
respectively. 

Let us close this subsection with the following observation. An interesting feature of the RG 
equation (|3.2ip is that it consists of a 'universal', /fc(-R)-independent part, and a truncation- 
dependent part. The former encompasses the terms S\,S k ,S k in the first line and captures the 
contribution of the scalars r], a, s, b, to and the auxiliary vector fields c M . The term S k thereby 
arises from combining the scalar traces with a different number of unphysical modes. The second 
part, which consists of the traces S k and S k , stems from the gravitational degrees of freedom 
K^ v and (j), and encodes all information about the particular form of fk(R)- Retrospectively, 
the appearance of this structure is not surprising: it is a direct consequence of the geometrical 
gauge-choice (|3.16p . which ensures that the gravitational degrees of freedom are carried by h^ v 
and 4> only. 

3.4 Heat-kernel techniques for evaluating operator traces 

The final step in constructing non-perturbative /^-functions from the flow equation (|3.2ip requires 
the explicit evaluation of the operator traces in terms of curvature invariants of the background 
manifold. This can be done by the standard heat-kernel techniques which are the subject of this 
subsection. 

In order to be able to apply heat-kernel methods to our problem, we first have to complete 
the traces by adding and subtracting the missing modes according to 



Here, W is an arbitrary, smooth, operator-valued function and s denotes the spin of the fields 
on which the covariant Laplacian acts. The discrete sum in (|3.24p can be evaluated using the 
multiplicities Di(d,s) and eigenvalues Ai(d,s) for —D 2 acting on a spin s field. To evaluate the 
complete traces, we assume that the function W(z), which arises from replacing the argument of 
the operator W(—D 2 ) by z, has a Fourier transform W(s). The trace is then expressed by this 
Fourier-transform, 





(3.24) 



le{h,-l m } 



Tr [W(-D 2 )] = lim / ds W(s) Tr e - {is ~ e) 

e ~*Q J — oo 



(3.25) 
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with e an infinitesimal regulator. The operator trace on the RHS can then be evaluated via the 
asymptotic heat-kernel expansion, 



Tr 



x d/2 

e -it(S 2 +Q) 1 ' ' 



\4TTt J 



J d d x^f§ [tra - ittia 2 (x,Q) - t 2 tia 4 (x,Q) + . . .] . (3.26) 



Here, Q denotes a "matrix potential" which is independent of the covariant Laplacian, the symbol 
tr denotes a trace with respect to the tensor indices of the fields on which the Laplacian acts, 
and the a2k(x,Q) are the heat-kernel coefficients of the Laplacian (see Appendix [A] for details). 
Substituting the heat-kernel expansion into (|3.25p . we obtain 

Tr [W{-D 2 )] ={^y d ' 2 j d d x^g {Q d /2[W}tva + Q d/ 2-i[W}tva 2 (x,Q) + ...} , (3.27) 
where the functionals Qnf^] are defined via 

/CO 
ds {-is + e)~ n W(s) . (3.28) 
-DO 

Making use of the Mellin transform, Q n [W] can be reexpressed through the original function 

W(z), 



Qn[W] = ' \ dzz n+l ~ 1 —^, »>-n,i€N. (3.29) 
r(n + i) J dz % 

For n > and n < G Z, it is convenient to choose i = and i = —n + 1, respectively. Then, 

assuming that lim z _ i . 00 W(z) = sufficiently fas1§|, (|3.29p simplifies to 

1 

Q n [W}=-— dzz n - l W{z), n>0, 

[n) J ° (3 30) 

d h W(z) 



Q_n[W]=(-l) 



dz n 



, h = —n G N . 

z=0 



This includes the special case Qo[W^] = W[0]. The second identity will play an important role 
when considering the full flow equation arising from (|3.2ip in Section [71 

For n > 0, the functionals Q n [W] resulting from (|3.21|) are closely related to the (generalized) 
threshold functions introduced in Appendix [Bj For the traces and <S|, this relation leads 

to the standard threshold functions (|B.ip , Defining 

N := (2g k y 1 d t (g k R k ) , (3.31) 

with g k an arbitrary function of the cutoff scale k, we have 

Qn[(Pk + c)- p Af] =k 2 ^ + V [$P( c /fe 2 ) + \d t \n{g k )¥ n {c/k 2 )\ . (3.32) 

Typically, g k represents a (combination of) /c-dependent coupling constants or wavefunction renor- 
malization factors multiplying R k . These enter the RHS of eq. (|3.32p through their "anomalous 
dimension" rj g (k) := dt ln(g k ). With regard to the later sections, it will be useful to introduce the 



3 For the functions appearing in the traces (|3.2ip . this holds. 
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anomalous dimension of Newton's constant and the non-local couplings v k and /j, k appearing in 
(g3D and (fL5|) . 

rj N (k) := -d t ln(Z Nk ) , r)y(k) := d t ln(v k ) and rjp,(k) := d t ln(/2 fe ) . (3.33) 

For the last trace in (|3.2ip . an analogous calculation relates Q n to the generalized threshold 
functions (|B.3|) . Using a schematic notation which highlights the — Independence of the trace 
argument, we find 



Qn 



8 t {g k (P 2 - D A ) + ~g k R k ) 

(u k P 2 + V k P k + W k )P 

= k 2 ^ 2 ^ [d t g k % 0>1 {u ky v k /k 2 ,w k /k A ) + Ag k T p nl (u k , v k /k 2 , w k /k 4 )} {3M) 
+ A ,2(n-2 P+1 ) ^dthTl^iu^Vk/k 2 ^,/^) + 2~g k r%{u k ,v k /k 2 ,w k /k 4 )} . 

With these tools, we have now all the ingredients to derive the /3-functions for particular choices 
of f k {R). To conclude this section, we illustrate the techniques just introduced by expanding 
the "universal terms" in the first line of (|3.21j) up to linear order in the curvature scalar. These 
results will be very useful when deriving the /3-functions of the ln(i?)- and i? _n -truncation in 
sections [5] and [H respectively. 

Let us first make a general observation on the "finite sum terms" in (|3.2ip . that is, on those 
terms which do not form a complete operator trace. For d > 2 and spherical backgrounds, 
these terms do not contribute to truncation subspaces which contain terms linear in R as the 
highest power of the curvature scalar. This is due to the fact that the trace terms are always 
proportional to the volume of the background sphere and, therefore, inherit extra inverse powers 
of the curvature scalar V oc Rr d l 2 . Thus, contributions coming from discrete eigenvalues are 
effectively suppressed by a factor R d l 2 and hence, for d > 2, only start to affect the running 
of couplings multiplying f d d x^/gR n ,n > 1. Given this, we can see that both 5| and those 
terms arising from the completion of traces in (|3.2ip do not contribute to the truncation spaces 
spanned by the ln(R)- and i? _n -truncation. Therefore, these can be ignored in the derivation of 
the corresponding /3-functions. 

With this observation, it is now straightforward to evaluate the contribution of the first line 
of the flow equation, neglecting 0(R 2 ) and higher curvature terms. Explicitly, 



S l = ~ (4^72 {°0 / d d xV9 + k- 2 ($-^% 2 + Cf ^ d/2 ^) J d d x^ , 
S ^ = ~ (i^72 {°o T / **y/9 + k- 2 (¥^/ 2 + Cf afy^) J d d x^R} . 



(3.35) 



The (i-dependent coefficients C| fc are defined in eq. (1A.12|) . In order to lighten our notation, we 
have introduced 

*£:=*£(<)), *£:=*£(()), T£. OT :=T£. m (l,0,0), := f £ ,(1, 0, 0) . (3.36) 
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Here and in the following, $n, <&n, Tn ;m and T^mi without argument will denote the cutoff- 
scheme dependent threshold constants defined above, while, for the related threshold functions, 
the argument will always be given explicitly. 

We close this subsection with the following remark. In order to evaluate the contributions 
of the traces f|3.25j) to our truncation, we use the early time expansion of the heat-kernel (|3.26p . 
expanding for small values of s. Based on perturbation theory, one would, at first sight, expect 
that an investigation of the gravitational RG flow in the IR should be based on the late-time 
expansion of the heat-kernel [57\ I58| . expanding the traces for large values s. Analyzing the 
Fourier-transformed functions W(s) arising from (|3,2ip . one however finds that, thanks to the IR 
cutoff, lim^oo W(s) = for k > 0\3 Thus, even for small values of k, the main contribution to 
the integral (|3.28|) arises from small values s, so that the early time expansion of the heat kernel 
can also be used to investigate the gravitational RG flow of the truncations (|1.4p and (|1.5p in the 
IR. 

4 General properties of the RG flow of /(i^)-gravity 

One of the most remarkable properties of the flow equation (|3.2ip is the fact that it is valid 
for any space-time dimension d, cutoff function R^, and arbitrary function fk(R)- This makes it 
worthwhile to pause and highlight some general features of the resulting RG flow, before analyzing 
concrete models specifying some of the quantities above. We start by formulating a decoupling 
theorem for non-local interactions in Subsection 14.11 and then move to derive the properties 
of fk(R) required for resolving the IR singularities as, e.g., the A = 1/2-line observed in the 
Einstein-Hilbert truncation, in Subsection 14. 21 

4.1 The perturbative decoupling of non-local interactions 

We now turn towards the discussion of the /3-functions for interaction monomials which are built 
from non-local curvature terms. Here, it is useful to distinguish two classes: interactions which 
blow up as the curvature scalar becomes small, as, e.g., R~ n - or ln(i?)-terms, and non-local 
interactions which remain finite as R — > 0, like i?ln(i?)-terms. Based on the flow equation (|3.2ip . 
one can argue that the /3-functions for the scale-dependent couplings multiplying interactions 
of the first kind are always trivial, while, in the second case, they are always proportional to a 
coupling constant multiplying an interaction monomial containing (some power of) m(i?), so that 
the corresponding interactions can be consistently decoupled from the RG flow. 

Let us start by investigating the /3-functions for the non-local couplings of the first kind. For 
this purpose, we will consider Tk[g] as a Laurent series in R, 

r*M = J dtxy/g Vn;kR n , n >0, (4.1) 

n>— no 

10 For the exponential cutoff function, W(s) even vanishes exponentially fast. 
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corresponding to fk(R) — WirGk X^n>-n n ^n;k 

Substituting this ansatz into (|3.2ip and 
dividing by V, we see that the LHS of the flow equation is a Laurent series whose highest order 
pole is given by R~ n ° . The coefficients in this expansion are proportional to the i-derivatives of 
the non-local couplings. 

Now, in order to extract the /3-functions for these couplings, we have to perform a Laurent 
series expansion of the RHS of the RG equation. Using the techniques introduced in Subsection 
13.41 this is done by first series expanding the arguments of the traces in R and then using the 
heat-kernel expansion (|3.26|) . ()3.27[) to evaluate the — D 2 -dependent traces. The first step reveals 
another remarkable property of the flow equation. Due to the homogeneity property of its RHS 
with respect to the function fk(R), expanding the arguments of the traces does not give rise to 
terms containing inverse powers of R. Thus, even though fk(R) contains poles as R — ► 0, the 
expansion of the trace arguments only yields terms which are regular as R — > 0. Taking into 
account that the heat-kernel expansion (|3.26p also contains positive powers of R only, we find 
that the RHS of the flow equation is regular in R = 0. Comparing the pole structure on both sides 
of the flow equation, we then conclude that the /3-functions of the non-local couplings appearing 
in (|4.ip are trivial, 

dtfin;k = , -n < n < , (4.2) 

verifying our first claim. Eq. (|4.2|) implies that non-local couplings of the first kind are scale- 
independent constants along an arbitrary RG trajectory. In particular, once we start with p, r = 
at the initial scale k, the RG flow does not dynamically generate non-local couplings of the form 
R~ n . 

The situation becomes slightly more complicated for non-local couplings of the second kind. 
In this case, expanding the arguments of the traces and using the heat-kernel expansion, one finds 
that such interactions also appear on the RHS of the flow equation. Thus, the /3-functions for 
this type of non-local couplings are non-trivial. However, since the heat-kernel expansion (|3.26l) 
does not contain (ln(i?)) ?1 -terms, the only source for the logarithmic factors in the non-local 
interactions is the expansion of the traces. This implies that the corresponding non-local terms 
on the RHS of the flow equation are necessarily proportional to a non-local coupling constant 
multiplying ln(i?)-terms in fk(R)- While it is difficult to construct the resulting /3-functions 
explicitly, we can still conclude that the resulting flow equation permits us to consistently set 
the corresponding non-local couplings to zero. In other words, if there are no interaction terms 
containing (powers of) ln(i?) present at the initial scale k, they will not be generated by the RG 
flow. Note that the same argument also applies to fractional powers of the curvature tensor, 
which also do not arise from the heat-kernel expansion (|3.26p , so that it is possible to consistently 
decouple such interaction monomials as well. 

This situation is reminiscent of the observations made in [13], where the RG flow of non- 
local functions fk(V) of the space-time volume V were studied. There, the /3-functions for the 
non-local couplings V 2 and Vln(V) were explicitly derived, with the result that, even though the 

As illustrated in the next section, this argument also applies to ln(i?)-terms which are not explicitly included 
in the ansatz (|4.1I) . 
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(3- functions of the non-local couplings were non-trivial, they could be consistently decoupled from 
the RG equations by setting these couplings to zero. Noting that such interaction monomials are 
also excluded in the heat-kernel expansion, it is tempting to speculate that all coupling constants 
whose interaction monomials do not appear in the heat-kernel expansion can be consistently set 
to zero in the flow equation. 



4.2 Resolving the IR singularities of RG trajectories with A > 

In this subsection, we discuss a feature of the RG trajectories found in the Einstein-Hilbert 
truncation (cfg. Figure [T]), namely, the termination of the type Ilia trajectories at Afc torm = 1/2 
for a finite value fcterm > 0. In principle, this class of trajectories would give rise to a positive 
cosmological constant in the IR, and it was argued in [56] that the RG trajectory realized by 
nature is of this type. In this light, it is an important question whether the termination of these 
trajectories is due to using an insufficient truncation to describe this part of the gravitational RG 
flow, and thus whether there are extensions of the Einstein-Hilbert truncation which allow the 
continuation of these trajectories down to k = 0. In the remainder of this section, we will use the 
RG equation for /(i?)-gravity to show that such extensions indeed exist. 

Let us begin by sketching the origin of the IR singularity in the Einstein-Hilbert truncation 



f k {R) 



-R + 2A fc . 



(4.3) 



From the RG equation (|3.2ip , we observe that its universal part is independent of A^ and therefore 
cannot be responsible for the singularity. Consequently, the singularity must arise from the f k (R)- 
dependent part. For conciseness, we illustrate its origin in noting that the analysis for 5| is 
completely analogous. 

Substituting the ansatz (|4.3p 



into S^, we find 



C.4EH 



Tr 



2T 



d t (Z Nk R k ) 



J Nk 



Pk 



(d-2)(d + l) kd 



• 2A k + 3(5=1) K 

<S> d/2 (-2A k /k 2 ) - ± m ¥ d/2 (-2A k /k 2 ) 



dt (ZNkRk) 
ZNk (Pk ~ 2A fc 



+ 0(R) 



(4.4) 



d d x^g~+0(R). 



2(4vr) d / 2 

Here, we carried out a Taylor expansion of the argument with respect to R in the first step 
and then used the heat-kernel expansion (|3.27|) . together with (|3.32[) . to obtain the last line. 
Substituting the threshold functions for the optimized cutoff (|B.9[) . this become J^l 

1 



4EH 



(d-2)(d+l) 
2(47r) d / 2 



r(d/2+l) 2T(d/2+2) 



k" 



1-2A* 



d?xy/g + <D(R) , 



(4.5) 



where we have expressed the cosmological constant through its dimensionless counterpart = 
A k /k 2 . From this expression, we see that the RHS of the flow equation, and therefore also the 
(3- functions for g k ,X k , has a pole at A = 1/2. An analysis of the RG flow in the Einstein-Hilbert 



12 For explicitness, we use the optimized cutoff from Appendix lB.2l Other cutoff-schemes like, e.g., the exponential 
cutoff or the sharp cutoff lead to qualitatively similar results. 
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truncation reveals that the type Ilia trajectories shown in Figure [I] run into the A = 1/2-line, 
which causes their termination at a finite value fcterm > 0. Our aim now is to characterize the 
properties of fk(R) which remove the poles at A = 1/2 from the RHS of the flow equation. Again, 
we will focus on St , noting that the same mechanism applies to 5|. 

The principal idea behind the resolution of the A = 1/2 singularity is to exploit the homo- 
geneity of the trace-argument in f k and its derivatives in such a way that the A^-terms arising 
in the denominators are multiplied by some power of R. The A^-terms are then dressed up as 
A k R e and thus expanded in the Taylor series. In that case, the RHS of the flow equation has 
no poles of the form (1 — 2\ k )~ p ,p > 0, and the singularity is resolved. This can be achieved if 
fl{R) contains a term which diverges in the limit R — » 0, 

lim f' k (R) = fl k R^ + ... , e>0, (4.6) 

H — >U 

with p, k a scale-dependent coupling constant. The relation (|4,6p is easily integrated to yield the 
corresponding condition on fk(R), 

]im f k (R) = R l - e + . . . , e>0, e + \ and lim f k (R) = v k ln(i?) + . . . , e = 1 . (4.7) 
R — >0 1 — 6 H — >0 

Let us illustrate the workings of this mechanism in detail by adding a term of the form (|4.7p 
to the Einstein-Hilbert truncation^ 



f k (R) = -R + 2A k + 167rG fc Jx k R 



l-e 



(4.8) 



Here, fi k has been rescaled for later convenience. Substituting this ansatz into S£, we find 



<7 4 



:Tr 2T 



z Nk d t 



^Nk \ — 



R e ) Ri 



^Tr 2T 



- R l +^ + 2k k Rt + 16vrG fc Jx k R 



*i/a(0) + i»7A*i/a(0)l / d d x^g + <D(R\R) . 



( 1J ^-m(Pk- 2 ^ 

_ (d-2)(d + l) d 
2(4vr) d / 2 

Using the threshold functions for the optimized cutoff, we observe that the poles at A = 1/2 have 
been removed. Thus, the singularities at A = 1/2 vanish with the inclusion of interaction terms 
of the form (|4.7j) . This will be further illustrated in the following sections, where it will be shown 
explicitly that the truncations (jl.4p and (|1.5p contain RG trajectories with positive cosmological 
constant in the IR Ao > 0. 

One cautious remark is now in order. Once terms of the form (14. 7p are included in the 
truncation subspace, one has to be very careful about taking the limit in which their corresponding 
couplings are sent to zero. The reason is that taking this limit does not commute with the Taylor 
expansion of the trace arguments. This can be easily seen from the following simplified example. 



13 For e = 1, one should replace R 1 t by ln(_R), which is easily done in the formulas below. 
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Setting g(x) = g\x 1 + g 2 ,h{x) = h\x 1 + h 2 , with gi,g2,hi,h 2 constant, and considering the 
expansion of the quotient 

g{x) _ g 1 x~ 1 + g 2 _ 9l + g 2 x _ 9l f g 2 9l h 2 \ 2 

h(x) ~ fcn-i + h 2 ~ ^ + h 2 x ~ h± l>i /i? ; 1 ; ' 1 j 

we see that taking the "non-local couplings" gi,hi to zero after the expansion leads to a quite 
different result than doing so before the Taylor series is constructed. From this example, it is 
clear that considering the RG equation on a truncation subspace which involves couplings of the 
form (|4.7p and setting the new coupling to zero does not automatically yield the /3-functions 
obtained in the Einstein-Hilbert case. That being said, in the remainder of this paper we will 
adopt the viewpoint that the /3-functions derived from the truncations (|l,4p and (jl.5p can be 
extended continuously to the point of vanishing non-local coupling. 



5 Truncation I: the In (R) -truncation 

Following up on our general discussion, we now continue our investigation by analyzing the non- 
perturbative /3-functions arising from the ln(i?)-truncation (jl.4p which, in the notation of Section 
El corresponds to 

f k (R) = -R + 2A k + 16irG k v k ln(R) . (5.1) 

We start by deriving the /^-functions governing the scale-dependence of the couplings G k ,A k 
and v k in Subsection 15.11 while their properties are investigated analytically and numerically in 
Subsections 15.21 and 15. 31 



5.1 The /3-functions of the ln(i?)-truncation 



Our first task in the construction of the /3-functions of the ln(i?)-truncation is the evaluation of 
the truncation-dependent traces S k and S k in (|3.2ip up to linear order in R. Substituting the 
ansatz (|5.1|) and using the heat-kernel techniques of Subsection 13.41 this yields 



<y4 

<->l. 



q 5 



(4,r)d/2 
k d - 2 

+ 



n 2T 
u 



( 47r )d/2 



{ 



a 



2T 



At 



' d-2 



d/2) 



+ c 2 2T 



*d/2-i + \ Vv *i/ 2 -i] } / dtx^gR , 



(4 
+ 



^i/2 { 2T d/2,i + h^d/2,0,1} j ddx V9 



kd ~ 2 f d + 2 ( 0T 2 

(47r)^/2l2(d-l) r ^ 
+ C 2 a/2-1,1 + hVv^d, 



k 2 (th? + Vn) j,i 
32irG k Vk 



d/2 



T 



d/2,0,0 + iVv (^d/2,1,1 ^d/2,0,0 



d/2-1,0,1 



d d x^g-R, 



(5.2) 
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for the 2T and scalar trace, respectively. Here, i]^(k) and r}y{k) are the anomalous dimensions 
(|3.33p and the C|, are defined in (|A.12p . Inserting the truncation ansatz into the LHS of the flow 
equation (|3.2ip in turn gives 

1 



d t T k [g] 



WttG; 



d a x^g -Rd t Z Nk + 2d t (Z Nk A k ) + / d*x^g ln(R) d t v k 



(5.3) 



Taking into account the contributions from the universal traces (|3.35p and equating the coefficients 
of the invariants J d d XyJ~g, f d d x^fgR and J d d x^fgln(R) spanning our truncation space then leads 
to the following system of coupled differential equations 



d t v k = , 
dtZNk 



W7rG k k d - 2 



(4 ^/2 _ 1^-1^1 cgT 



+ 



- d 2 -6 
, 6d 



d/2 

1 

'd/2- 



16nG k v k 
2 



r^W d+2( T 

°2 )^d/2-l ~ d-l\ L d/2,2 



d(d-l) 



d/2 



(d-2)(d+l) A fc _ d 3 -4d 2 +d-3 \ $2 

XT 4 ! ^ XT 1 ^ 

2 L d/2,l) 3 x d/2-l,l 



(5.4) 



f,//Cd T [(rf 2 -3d-2)^ /2 + 4Ti /2il ] 



(4vr) d /2- 

The /3-functions resulting from these equations are most conveniently expressed in terms of the 
the dimensionless couplings 

g k = k d ~ 2 Z„\ G k = k d ~ 2 G k , A fc = k~ 2 A k , v k = k~ d v k 



(5.5) 



In terms of these, the eqs. (j5.4j) imply the following set of autonomous first order differential 
equations 

9t9k = @g(g,\,v), d t X k = Px(g, X,v) , d t v k = v (g, X,v) , 



(5.6) 



where 

Pg(9, X,v) = (d-2 + t)n) 9k , 

Px(jg, X, v) = - (2 - VN ) Afc + (4vr) 1 - d / 2 5fc ((d 2 - 3d - 2) <S> l d/2 + 4 Tj /2>1 
0vig,X,v) = -dv k . 
Here, the anomalous dimension of Newton's constant is given by 



(5.7) 



VN 



32irg k 



2(A^) d / 2 -v- k ^ d/2 C 2 ^ 



- (rf-2)(rf+l) A fc _ d 3 -4:d 2 +d~3 \ fo2 
. Wng k v k d(d-l) ) ^d/2 



(5.8) 



+ 



' d?-% 
, 6d 



C 2 > 



rf+2 ^2 l-y-l \ 

I Mfti o L d/2,l) 



d/2— 1 d-lV J 'rf/2,2 2" L d/2, 



3 1 d/2-l,l 



This completes the derivation of the /3-functions for the ln(Z2)-truncation. Note that these (3- 
functions are consistent in the sense that they incorporate all contributions proportional to the 
interaction monomials spanning the truncation subspace. As expected from the general discussion 
in Section I4TT1 the /3-function for the dimensionful non-local coupling v k is indeed trivial, dtv k = 0. 
Furthermore, the only poles of (|5,7p occur at the line 



I(4vr)-^C 2T ^ 



d/2' 



(5.9) 
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showing explicitly that the singularity at A = 1/2 has been resolved, in full agreement with the 
arguments given in Section HI 

A second property which can be directly deduced from the /3-functions (|5.7p is the decompo- 
sition of the coupling constant space in various sectors whose trajectories do not mix under the 
RG flow. Since j3 g {g = 0, X,v) =0 and @ v (g, A, v = 0) = 0, the g = 0- and v = 0-planes cannot 
be crossed by the trajectories. In other words, a RG trajectory starting with positive g^ (or vt) 
will never evolve to gj. < (or Vk < 0), and vice versa. 

5.2 The fixed points of the RG flow 

In order to gain some first insights on the dynamics arising from (|5.6p . it is very helpful to 
analyze the fixed point structure of the /3-functions. At the fixed points (FP), all /3-functions 
vanish simultaneously, 

P g (g*,X*,v*) = 0, /3 A ( 5 *,A>*) =0, (3 v (g* , \* ,v*) = 0, (5.10) 

and many properties of the RG flow can be captured by linearizing the RG equation at such points 
and studying the flow in their vicinity. Introducing the generalized couplings Uk = [fffc, Afc, t^] 
and the Jacobi matrix 

RR-. 

(5.11) 



B = {Bij) , B^ - 



duj 



FP 



the linearized flow equation reads 

dm « ^Bijtiij - u*) . (5.12) 

3 

Using the stability coefficients 8 l = —A 4 , with A* the eigenvalue of B associated with the right 
eigenvectors V 1 , B^ can be diagonalized, so that eq. (|5.12p is easily solved. These analytic 
solutions then provide a good picture of which RG trajectories are dragged into the FP or repelled 
along an unstable direction. 

Applying this construction to the /3-functions (|5.7p . we first observe that j3 v = forces all 
FP onto the Vk = plane. Substituting this condition into the remaining /3-functions, (5 g and (3\ 
simplify considerably and their roots are easily found analytically, 

GFP : {g* = 0, A* = 0, v* = 0} , 

IRFP: {g* = 0, \* = ^-,v* = 0}, (5.13) 

r * (47r) d / 2 - 1 d(d-2) 5* , „ (rf-2) 0* , 

NGFP: {g* = -j- r , A = — - .1 d/2 , v* = } . 

Here, GFP denotes the Gaussian Fixed Point, while the other FP have providently been labeled 
Infrared Fixed Point (IRFP), and Non-Gaussian Fixed Point (NGFP). Let us then discuss the 
stability properties of these fixed points in turn. 
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At the GFP, the stability matrix B becomes 

d-2 



B 



GFP 



(4vT 



,i-d/2, 



; ((d 2 -3d-2)d>i /2 + 4Ti /21 ) 


with stability coefficients and associated right eigenvectors 

V 1 




-2 
-d 



(5.14) 



2 
d. 



d, 



l,i(4 7 r) 1 -^((d 2 -3d-2)^ /2 + 4Ti /21 ),0 



^ 2 = [0,1,0] T , 

v 3 = [o,o, 1] 



"d/2, 



(5.15) 



T 



This information can be used to solve the RG equation in the vicinity of the GFP, giving 

k d - 2 



9k =«i 

Afc =OL2 

Vk =a 3 



M d~2 
k- 2 



k 



d-2 



M- 

k- d 



(A \\-d/2 

2 + «i -^j— ((* - 3d - 2)*^ + 4T i /2jl ) ^ 



(5.16) 



Here, {ctj} denotes the set of integration constants and M is an arbitrary fixed mass scale. It is 
then natural to definite the dimensionful integration constants 



Go = ml\ d , 



A = a 2 m P1 , v = a 3 mpj , 



(5.17) 



where we have identified M with the Planck mass mp\ by setting a\ = 1. By using the relations 
()5.5p . we then obtain the scale-dependence of the dimensionful coupling constants in the vicinity 
of the GFP, 

Gk = Gq , 



Ai 



(A \l-d/2 

A + L-L Go ( (d 2 _ M _ 2 )$i /2 + 4Tj /3)1 ) 



(5.18) 



v k =v . 

These equations show that Gfc,Afc and Vk run towards constant, non-zero values as k — * 0. In 
particular, A^ does not vanish in the IR unless we make the special choice Q2 = 0. Thus the GFP 
does not determine the behavior of the cosmological constant Ao in the infrared. 
Let us now turn to the IRFP (|5.13p . Here, the matrix B becomes 



B 



IRFP 



d 







(4^)^/ 2 ((d 2 -3d-2)^ /2 + 4Ti J 2 



d/2 ^ d/2, 







4(47r) d / 2 
(d+l)(d-2)$2 /2 

-d 



(5.19) 



and the eigensystem of Birfp is 



-d, 
-2, 



r 1 



1 , ^ (4^-V 2 ((d 2 - 3d - 2)*J /2 + 4Ti /2a ) , 



"d/2, 



^ 2 = [0,1,0] T , 



0. 



d/2 



4(4tt) 

(d 3 +d 2 -4d-4)$jj /2 ' 



1 T 



(5.20) 
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Note that here the stability coefficient #2 has now changed sign compared to the eigensystem of 
the GFP (|5.15p . making the eigendirection V 2 IR-attractive (hence the name IRFP). Solving the 
linearized flow equation results in 

k d 

p (a \l-d/2 t.d 

X k =*2 W2 + ax LJ_ -3d- 2)^ /2 + 4T$ A1 ) w 

4(47r) d/2 k -d $l /2 ( 5 - 21 ) 

" 3 (d 3 + d 2 - 4d - 4) $ 2 /2 M" d + 2$ 2 /2 ' 

again denoting the integration constants by on. Using the relations (|5.5|) . we then obtain the 
scale-dependence of the dimensionful coupling constants, 

2 



k 

Gk = Go — t - ' 



u4 (a \l-d/2 h.d+2 

A fe = Ao-^ + ((d 2 - 3d - 2)^ /2 + 4TJ A1 ) Gb-3- 



PI " - '"PI 

2 



(5.22) 



Thus, for d > 2, the IRFP has an unstable direction in driving the RG trajectories with 
i)o/0 away from the FP. For trajectories in the plane vq = 0, however, the IRFP is an infrared 
attractor for both Newton's constant and the cosmological constant. In this case both Gk and 
Afc are dynamically driven to zero as k —* 0, independently of their values at the initial scale k. 

Finally the NGFP is located in the (probably unphysical) region g < 0. Determining its 
stability coefficients, one finds that it is UV attractive for all three couplings gk,^ k and vp.. In 
this sense, the properties of this NGFP are very similar to the ones observed in the Einstein- 
Hilbert truncation, cfg. Fig. [Q except that its position has now shifted from g* > 0, A* > to 
(/* < 0, A* < 0. We will return to the discussion of this particular feature in our concluding 
remarks in Section [HJ 

5.3 Numerical solutions of the flow equations 

After discussing their fixed point structure, we now proceed by integrating the flow equations 
(|5.6|) numerically. This requires fixing the space-time dimension and the cutoff-scheme explicitly, 
so we will restrict ourselves to d = 4 and to the optimized cutoff (|B.7p in the following. Given the 
decomposition of the coupling constant space, we will furthermore limit our focus on trajectories 
with gk > 0, as these are the physically most relevant ones. For d = 4 and the optimized cutoff, 
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Figure 2: RG flow of the ln(J?)-truncation in the fixed plane v = 0. The flow is dominated by the GFP at 
g* = 0, A* = and the IRFP at g* — 0, A* = 1/6. All RG trajectories running to the left of the "separatrix" are 
captured by the IRFP as k — > 0. 

the FP ([5TT3]) are located at 

GFP : {g* = , A* = , v * = 0} , IRFP : {g* = , A* = ± , v* = 0} , 

(5.23) 

NGFP : {g* = -?f , \* = -1 , v * = 0} . 

We note, in particular, that the NGFP is indeed situated in the region g < and hence will not 
play a role in the subsequent discussion. Moreover, the singular line (|5.9p is located at 

v S i ng = 0.00263 . (5.24) 

Due to the existence of the IRFP, the most interesting region of the parameter space is the 
one comprising the RG flow of g^, \k ° n the fixed plane v = 0, which is shown in Figure (2) This 
diagram clearly illustrates the separation between trajectories with positive and negative values 
gt- The RG flow in the upper half plane is thereby dominated by the GFP and the IRFP. The 
RG trajectories in this region fall into three distinct classes, which, in analogy to the Einstein- 
Hilbert truncation shown in Figure [lj are labelled Type la, Ha, and Ilia, respectively. The class 
Ha consists of the single trajectory ("separatrix") hitting the GFP as k — > 0. The trajectories 
of Type la run to the left of the separatrix and lead to a negative cosmological constant Ao in 
the infrared. The trajectories of Type Ilia, to the right of the separatrix, are captured by the 
IRFP as k — > 0, which implies that the corresponding dimensionful couplings Gfc,Afc go to zero 
as k — > 0. 
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The IR behavior of these classes can also be understood on the level of the linearized solutions 
(|5.18p and (|5.22p . as follows. The classes la and Ha correspond to linearized solutions «2 < 
and Q2 = 0, respectively. The solutions with «2 > are driven away from the GFP regime and 
captured by the IRFP f|5.22|) . giving Go = 0, Ao = independently of any integration constants. 
The latter feature is also illustrated in Figure [3] (a). Due to the absence of an UV FP in the ln(i?)- 
truncation, we do not expect, however, that the classes of RG trajectories discussed here give 
rise to a well-defined behavior as k — > oo. This is possibly owed to the fact that our truncation 
subspace is insufficient for a proper description of quantum gravity in the UV. 

Going away from the v = 0-plane, we have to distinguish the two instances v k < and 
v k > 0. To illustrate the behavior in these regions, some examples for trajectories starting with 
A(k = 1) > are shown in Figure El For the trajectories with v/- < we find that both G k ,A k 
decrease with decreasing values k. Comparing the series of diagrams (a), (b), (c), in which the 
modulus of the non-local coupling is continuously increased, we observe that this "quenching" 
works more efficiently for small non-local coupling \v k \, the case of v k = 0, corresponding to the 
IRFP, being the most efficient. 

Finally, for the trajectories with v k > 0, we distinguish two different behaviors. Trajectories 
starting with vt < f s ing run into the singularity and terminate at finite fcterm = (t^/using) &j 
whereas trajectories in the region vt > losing are non-singular and can lead to increasing values of 
Gk, Afc as k decreases. An example of this latter behavior is shown in Figure (d). 

6 Truncation II: the ^""-truncation 

We now turn towards our second class of non-local truncations, the i? _n -truncations (|1.5|) . for 
which 

f k {R) = -R + 2A k + 16irG k fL k R- n . (6.1) 

We begin by deriving the non-perturbative /3-functions of this truncation in Subsection 16.11 and 
subsequently discuss their properties in Subsections 16.21 and 16.31 

6.1 Deriving the RG equations 

Following the calculation of the previous section, we first compute the contribution from the 
truncation-dependent traces Sf, and 5|. Substituting the ansatz (|6.ip into (|3.2ip and evaluating 



32 




t.AU ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' K U.Z.W ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



(c) (d) 

Figure 3: Numerical solutions to the ln(7?)-truncation flow equations in terms of the dimensionful couplings, for 
(a) Vj, — 0, (b) fjj = -5 x 1CF 5 , (c) v k = — 5 x 1CP 2 and (d) v k = 0.01 imposed at k = 1. Here, we have taken 
G k = 0.25, while A k assumes various values. We see that, for v k < 0, the cosmological constant is quenched in the 
IR limit, going to zero (a), or small but non-zero (b) values. 
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the resulting expressions up to linear order in R, we find 
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Here, r]n{k) is the anomalous dimension (|3.33p . and the cutoff-scheme dependent constants are 
defined in (|3.36p . Conversely, the LHS of (j3.21[) now reads 

d t f k [g] = — l — I d d x^[ - R8 t Z Nk + 2 d t (Z Nk A k )\ + I d d x^^ d t fl k . (6.3) 
k J J 

Comparing the coefficients multiplying the interaction monomials in the truncation ansatz again 
leads to three coupled differential equations encoding the scale-dependence of Zj^ k ,A k and ji k . 
Rewriting these equations in terms of the dimensionless coupling constants 



g k = k d ~ 2 G k , X k = k' 2 A k , n k = k-( d+2n ^ji k , 
we then arrive at the RG equations of the i?~ n -truncation, 

9t9k = Pg{g,K m) > d t Xk = P\{g, A,/i) , d t fj, k = ^(g,X,ti) , 



(6.4) 



(6.5) 



with 



Pg(g, A, y) = (d - 2 + ri N ) g k 



x (g, A, M ) = - (2 - m) Afc + (4vr) 1 - d / 2 5fc ((d 2 - 3d - 2) *J /2 + 4 T l d/2 l 
Pn(g, A, /i) = - (d + 2n) M fc . 
The anomalous dimension of Newton's constant is given by 
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The d-dependent constants C| fe are given in (|A.12j) . Observe that H nj d is actually a (truncation- 
dependent) constant, independent of any coupling constants, so that r/jv is linear in g^. This 
feature manifestly differs from the rjjy obtained in the Einstein-Hilbert truncation [3J, which con- 
tained contributions from all orders in g^. Thus, the i? _n -truncations probably do not capture 
non-perturbative contributions from an infinite number of graviton loops. Moreover, again con- 
firming the results of Section 0J we note that the /3- function for the non-local coupling fik is trivial, 
so that fik is a fe-independent constant along a RG trajectory. We thereby notice that the non- 
local coupling /Ufc actually decouples from the flow of gk, A&, in the sense that g (g, A, (j,),(3 g (g, A, Li) 
are independent of Hk- Thus, the RG flow for gk,^k is not altered by the value of the constant 
fik and can be studied independently from the flow of ii^. 

Finally, we observe that the /3-functions have no poles in the coupling space, being well-defined 
for any value g, A, li. Similarly to the ln(i?)-truncation, the coupling constant space decomposes 
into various regions, whose RG trajectories do not mix under the RG flow. These regions are 
separated by the fixed planes g = and \i = 0, as well as the fixed plane at g* = (d — 2)/H n ^. 
The latter is due to the linearity of r]N, which implies that f3 g (g*, A, li) = for any value A, fj,, so 
that the g*-plane cannot be crossed by the RG flow either. 

6.2 Fixed points of the RG flow 

In order to get a better idea of the RG flow encoded by (|6.6|) . we follow the strategy of Subsection 
15.21 and analyze the FP of the /3-functions. The conditions (3u(g* , A* , /x*) = implies that all FP 
must be located on the li = 0-plane. Solving f3 g (g*, A*, fx*) = 0, (3\(g*, A*, /x*) = analytically, we 
find that there are two roots, yielding a GFP and a NGFP, 

GFP : {g* = 0, A* = 0, Li* = 0}, 
NGFP: { g * = -pi , A* = -^— [ - J-*£ ™,^* = 0}. (6.9) 

Focusing on the NGFP for the moment, we see that its position depends crucially on the sign of 
H n ^ which, for fixed d, in turn depends on the exponent n of the i? _n -term in our ansatz. Noting 
that g* < (g* > 0) for H n ( i > 0, (H n ,d < 0), respectively, we define n cr i t as the (possibly non- 
integer) positive real root of the quadratic equation H n .d = 0, assuming it exists. The existence 
and actual value of the transition point n cr i t are thereby cutoff-scheme dependent and hence not 
universal features of our flow equations. Using the optimized cutoff, this particular feature of the 
i?~ n -truncation will be further investigated in the next subsection. 

Let us now proceed by discussing the stability properties of the GFP and NGFP in turn. 
Linearizing the RG flow at the GFP, the stability matrix (|5.1ip becomes 



B 



GFP 
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(4^)^ 3d -2)^ + 4^) -2 

-d-2n 



(6.10) 



35 



with corresponding stability coefficients and respective right eigenvectors 
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M-2)*J /a + 4Tj /J , |1 ),0 



(6.11) 



With this information at hand, it is then straightforward to solve the linearized flow equations 
close to the GFP. The result is here very similar to our findings in the ln(i?)-truncation, indicating 
that Gk, Afc and fik approach a^-dependent constants as k — > 0. In particular, the IR value of 
vanishes only if the corresponding integration constant c*2 = 0. 

Turning to the NGFP, we find that the stability matrix B becomes 
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Looking at the stability coefficients 0j, we note that, for d > 2, the NGFP is UV attractive along 
all eigendirections V % . Solving the linearized flow equation, we then find 

k 2 - d d-2 



9k =OL\ 



M 2 ~ d H,, 
k- d (4vr 
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l-d/2 



' M~ d 



d 2 



((d 2 - 3d - 2)& d/2 + AT') 



k 



2-d 
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(6.14) 



M~ 
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where M again denotes a fixed mass scale. This result indicates that the NGFP is indeed a UV 
attractor for the RG flow, i.e., for k — > oo, all the dimensionless coupling constants gk^k^k 
approach their fixed point values (|6.9[) independently of the integration constants «j. In the case 
of H nt d < g* >0,A* >0 this FP is the natural generalization of the NGFP shown in Figure 

m 



6.3 Numerical solutions of the flow equations 

After discussing the general FP structure of the /3-functions (|6.6p . we now restrict ourselves to 
d = 4 and the optimized cutoff (|B.7p and proceed with the numerical investigation of the RG 
flow. In this case, the FP (|6.9p are located at 

GFP : {g* = , A* = , // = 0} , NGFP : {<?* = ^ , A* = ^ , // = 0} , (6.15) 
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(a) (b) 

Figure 4: Coordinates of the NGFP on the fit = 0-plane for different values of n, using d — 4 and the optimized 
cutoff. The transition from #ngfp < to (?ngfp > occurs for n C iit — 9.09. 



with the function (|6.8p determining the position of the NGFP, 

_ 7n 2 + 57n + 60 
= 24vrn(n + 1) " (6 - 16) 

We observe that H ni 4 as a function of n, indeed has a zero at positive values n cr it — 9.09. At this 
value the NGFP passes from < 0, A* < at n < n crit to > 0, A* > for n > n cr i t . This 
n-dependence of the NGFP is illustrated in Figure [H Thus, for n > n cr it, the i? -n -truncation 
indeed gives rise to a NGFP which is a UV attractor for g^ , A^ , /i^ . Based on this observation, 
we will now discuss the RG flow for truncations with n < re cr it and n > re cr it in turn. 

Let us first consider the case n < n cvlt . Choosing n = 1, the typical RG flow in the [i = 0-plane 
is exemplified in Figure [5j This diagram nicely illustrates the separation between trajectories with 
negative and positive g^. The RG flow on the upper half plane, g > 0, is completely dominated 
by the GFP. In analogy with the previous section, we can again classify the RG trajectories in 
this region as type la or Ilia, according to whether they lie to the left or to the right of the 
"separatrix" , i.e., the RG trajectory of type Ila hitting the GFP as k — > 0. We note that the 
trajectories of all three classes can be continued to the IR, k — ► 0, where they give rise to a 
negative, positive or vanishing value of respectively. Due to the absence of a UV FP, it 
is unclear, however, whether these trajectories give rise to a well-defined theory in the UV as 
k — > oo. 

For n > n cr j t the RG flow in the upper half plane is governed by the interplay of the NGFP 
and the GFP. Using the example n = 100, this is illustrated in Figure Here, we see that all 
RG trajectories with gj. > emanate from the NGFP in the UV. They can then be classified into 
A-type and B-type trajectories, depending on whether they run below or above the ^ GFP -plane. 
According to their behavior in the IR, these trajectories can furthermore be divided into type la, 
Ila and Ilia, depending on whether they give rise to a negative, zero ("separatrix") or positive 
cosmological constant Aq. The remarkable feature of this phase diagram is that all trajectories 
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Figure 5: Typical RG flow of the R "-truncation for n < n cr it on the fixed plane fik = 0, choosing n — 1 for 
explicitness. 

with gk > have a well defined IR limit k = 0. Taking into account the //^-independence of the 
flow, all RG trajectories on the upper half plane are complete, in the sense that they give rise to 
a well-defined Tk[g] on the entire interval < k < oo. 

The typical RG flow of the dimensionful couplings Gfc,Afc along trajectories of Type Ilia 
(n < n cr ;t) and A- Type Ilia (n > n cr ; t ) is then shown in the diagrams (a) and (b) of Figure 
respectively. In both cases we find essentially no significant decrease of the cosmological constant 
as k — > 0. Moreover, due to the sign difference in its /3-function we find that Gk increases 
(decreases) towards the IR. The RG flow of B-type trajectories is qualitatively very similar to the 
one found in the A-region and therefore not shown explicitly. 

7 A differential equation capturing the RG flow of fk(R) 

This section adopts an orthogonal approach to the truncations previously discussed, in the sense 
that we are not restricting ourselves to a specific form of the function fj-(R), but instead specialize 
to d = 4 and a special choice of cutoff-scheme, the optimized cutoff (|B.7p . Using the heat-kernel 
results collected in Appendix [A] then allows us to evaluate the traces in (|3.2ip explicitly for 
arbitrary function fk(R), resulting in an autonomous partial differential equation governing the 
RG flow of /(-R)-gravity. 
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L 10 20 30 

Figure 6: Typical RG flow of the i? _n -truncations with n > n C iit, illustrated by the choice n = 100, in the fixed 
plane jj-k ~ 0. The flow is dominated by the interplay between the NGFP at g* > 0, A* > and the GFP at the 
origin. All RG trajectories are well-defined at all scales < k < do. 

7.1 Deriving the non-perturbative flow equation 

We start with the following observation F 4 ] When the operator traces appearing in (|3,2ip are 
evaluated via the heat-kernel expansion, their arguments W(z) enter into the series expansion 
via the functionals Q n [^]- In d = 4, the coefficients multiplying the heat-kernel coefficients 
tr[ti2i], I > 2, are thereby proportional to Q2-z[I^] which are essentially the Ith derivative of W(z) 
evaluated at z = (|3.30[) . The virtue of working with the optimized cutoff appears when (|B.7p is 
substituted into W{z) and the derivatives are computed explicitly. Then, one finds that starting 
from a finite value /term ah derivatives vanish at z = Oj^jj Thus, working with the optimized 
cutoff has the major advantage that the knowledge of a finite number of heat-kernel coefficients 
suffices to evaluate the traces (j3.21|) . We stress that we do not expect this feature to persist once 
a different cutoff-scheme as, e.g., the sharp cutoff or the exponential cutoff, is used. 

For the purpose of this section, it is useful to redefine fk(R) by absorbing the wave function 
renormalization Z^k, working with fk(R) '■= Zj^kfk(R) and dropping the tilde in the following. 
Furthermore V denotes the volume of S 4 . 

Equipped with the result above, we now proceed with the explicit evaluation of the RHS 
in (|3,21|) , In contrast to the truncations discussed previously, we now do not perform a series 

14 We are in debt to R. Percacci for explaining this point to us. 

15 Strictly speaking, this "derivative expansion" truncates for the traces Si and <S| only. The traces and 5^ 
give rise to an infinite series of terms which, starting at itorm are proportional to 8{k 2 ) and its derivatives. In the 
following we will neglect the contributions of such terms to the flow equation. 
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Figure 7: Numerical solutions to the flow equations of the i? _n -truncations for (a) n < n cr it and (b) A-type 
trajectories with n > n ct it, in terms of the dimensionful couplings. Imposing the initial conditions at k — 1 we 
have set Gj, = 0.25 with At taking various values. Due to its decoupling these results are valid for any value 



expansion of the trace arguments with respect to R. Evaluating the terms in the first line then 
gives 
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Here, the terms proportional to V arise from the evaluation of the complete traces, the heat-kernel 
coefficients a2k\s being given in (|A. 10|) and (|A.11|) . Furthermore, the <&n(w) denote the threshold 
functions for the optimized cutoff (IB.9h . The V- independent terms capture the contribution of 
the finite number of — D 2 -eigenmodes required to complete the traces (cf. Subsection [32 



The evaluation of the /fc(i?)-dependent traces S k and <S| is slightly more involved. For the 
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2T trace, we find 
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with heat kernel coefficients a 2 k\2T listed in (jA.lip . The argument of the threshold functions is 
given by the dimensionless quantity 
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(7.3) 



and, in analogy to (|3.33p . we have defined the "anomalous dimension" of f' k (R) as 

Vf (k):=d t Hf' k (R)). (7.4) 

Finally, we evaluate the scalar trace S k . Setting d = 4 and comparing its argument W$(z) 
with the schematic notation used in (|3.34|) motivates the definition 

u k := mf k ' , v k := 12k~ 2 (f' k - 2Rf' k P ) , w k := 4k~ 4 (2f k - 2Rf k + R 2 f, 

In terms of these, the derivatives of W5 (2)^=0 become 
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As mentioned above, the higher order derivatives Wg"'(0) are proportional to 5(k 2 ) or derivatives 
thereof and will not be included in the following. Using the generalized threshold functions (|B.3|) . 
we can then explicitly evaluate S k , 
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with the heat-kernel coefficients a2fc|o from ()A.10p . We thus see that, thanks to the optimized 
cutoff, the RHS of eq. (|3.2ip can indeed be explicitly computed for general function f k (R). 
Substituting the /fc(i?)-ansatz into the LHS, on the other hand, yields 



d t Tk[g]=2K 2 Vd t fk(R). 
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The flow equation governing the scale-dependence of f k (R) is then obtained by substituting 
dUD , (1772]) and ([777]) into the RHS of fl372"T|) and equating this to £L8|). Explicitly substituting the 
(generalized) threshold functions from Appendix lB.2l then leads to the following partial differential 
equation for f k (R), 
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In order to discuss its properties, it is useful to rewrite eq. (|7,9|) in terms of dimensionless 
quantities. This is done by first noting that the volume and curvature scalar of the background 
4-spheres are related by V = 3847r 2 i?~ 2 . Using this relation to eliminate V, eq. (|7.9p becomes 
a partial differential equation in k, R. The dimensionful quantities R and fk(R) are then traded 
for their dimensionless counterparts, 



p := k 2 R and ^Fk(p) '■- 



1 



16vrGi 



k- A f k (R/k 2 



In terms of these, 



dtfk = 16vrG^ k 4 (d t T k + 4F k - 2pT' k ) , 
while the anomalous dimension (|7.4p becomes 
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and similar relations hold for the other derivatives of f k (R). Reexpressing the dimensionful 
quantities in (|7,9p through the dimensionless ones leads to the following autonomous partial 
differential equation governing the RG flow of J~ k (p), 
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This equation then constitutes the desired autonomous partial differential equation for the 
RG flow of /(H)-gravity. 

A very peculiar feature of (|7.13p is the appearance of discontinuities on its RHS, induced 
by the 0-function terms. The origin of these contributions can be traced back to the use of the 
optimized cutoff, for which the finite sums of the type (|3.24p give rise to step functions. The 
observed discontinuities then probably reflect the price to be paid for being able to truncate the 
heat-kernel expansion at a finite order. Repeating the construction above using a different cutoff- 
scheme will eventually remove these stepfunctions but requires one to carry out the heat-kernel 
expansion to all orders in R. We also note that similar discontinuities have already been observed 
when constructing the /3-functions of the i? 2 -truncation using the exponential cutoff-scheme [9] . 
In that case, the discontinuities were not in the parameter p but in the space-time dimension d. 
Evaluating the finite sums for the exponential cutoff induced contributions to the flow equation 
which were proportional to and 5^d, so that the flow equations could change discontinuously 
with the space-time dimension. While it is clearly desirable to get a better understanding of such 
discontinuities, we will nevertheless leave this topic for the time being, and rather elucidate some 
other properties of the flow equation (|7.13p . in the next subsection. 

7.2 Fixed points of the RG equation 



From the dimensionless flow equation (|7.13p . it is now straightforward to obtain an ordinary 
differential equation describing the fixed functionals J~*(p) of the RG flow of /(-R)-gravity. By 
definition, these satisfy dt^kip) = and are consequently solutions of the non-linear ordinary 
differential equation 
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Due to its non-linearity, finding exact analytic solutions to this equation is rather involved, and a 
detailed analysis of the RG flows following from the partial differential equations (|T. 13|) and (|7.14p 
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0.00523 


-0.0202 
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n nnQQQ 
U.UUooo 


-U.Ulzo 


U.UU14y 
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0.00518 


-0.0196 


0.00070 


-0.0104 








4 


0.00505 


-0.0206 


0.00026 


-0.0120 


-0.0101 






5 


0.00506 


-0.0206 


0.00023 


-0.0105 


-0.0096 


-0.00455 




6 


0.00504 


-0.0208 


0.00012 


-0.0110 


-0.0109 


-0.00473 


0.00238 



Table 2: Location of the NGFP obtained from the fixed functional equation (|7.14p by expanding 
Tk{p) in a power series in p up to order p n , including ^-dependent coupling constants Ui(k), 
i = 0, . . . , n. 

is beyond the scope of the present paper. Nonetheless, let us end our discussion by pointing out 
some of their properties. 

We first investigate the possibility of a GFP arising from the flow equation (|7.14p . In this 
course, we reinstall the dimensionless Newton's constant by setting 

jr *(p) = 7K—F*(p)> ( 7 - 15 ) 

where g* indicates that g^ is taken at the fixed point, dtgu = 0. Substituting (|7.15p into (|7.14j) . 
we see that, due to its homogeneity properties, the RHS of the flow equation is independent of 
g*. Schematically, (|7.14p then takes the form 

768vr 2 - p T'\ = 16tt 5 * [...], (7.16) 

with the terms inside the bracket being independent of g*. At the GFP, we have g* = 0, so that 
the "quantum corrections" on the RHS decouple, and (|7.16p reduces to 

2^ GFP -p^ GFP = 0. (7.17) 

This equation has the solution ^p FP = (3*p 2 , with (3* an arbitrary integration constant. The 
corresponding fixed functional is given by 

f GFP = T7T~ f dtxy/ghie . (7.18) 

167TSI* J 

We note that, with respect to classical power counting, this is just the marginal operator of 
/(-R)-gravity. Furthermore, for (3* finite, the coupling /3* = (lGirg*)^ 1 (3* gets shifted to infinity 
as 3* — ► 0, which is in agreement with the vanishing of the GFP observed in [9]. 

Following [22], the flow equation (|7.13p and (|7.14p can also be used to investigate the UV 
properties of the RG flow. This has a close relation to the asymptotic safety conjecture of 
quantum gravity, as it allows one to (perturbatively) test the existence and stability properties of 
the NGFP in truncation spaces incorporating successively higher powers of the scalar curvature. 
In this light, we consider the UV limit of (|7.13p . For fixed background curvature, this corresponds 



44 



n 


9' 


nil 

9 


02 


03 


#4 


05 


06 


1 


2.38 


-2.17 












z 


l.ZO 


O A A 


z ( .U 










3 


2.67 


-2.26 


2.07 


-4.42 
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-2.42 
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-4.28 


-5.09 
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2.57 


-2.67 


1.73 


-4.40 


-3.97 + 4.57i 


-3.97 - 4.57i 




6 


2.39 


-2.38 


1.51 


-4.16 


-4.67 + 6.08i 


-4.67 - 6.08i 


-8.67 



Table 3: Stability coefficients for the NGFP with increasing the dimension of the truncation space 
n + 1. The first two critical exponents are a complex pair 9 = 9' + i9". 

to expanding (|7.13[) for p = R/k 2 <C 1. In this limit, all ^-functions have positive arguments and 
thus contribute to the flow equation. 

We then make a polynomial ansatz for J~k{p), setting 

n 

Fk{p) =Y, Ui( V P l ' re >!eN. (7.19) 

i=0 

Here, u\{k) = — (167r<?fc) , so that a negative u\(k) actually corresponds to a positive Newton's 
constant. Substituting this ansatz into (|7.13|) . expanding its RHS in powers of p, and matching 
the coefficients up to p n yields the non-perturbative /3-functions for the dimensionless couplings 

Ui(t), 

d t Ui = (3 Ui (u , . . . ,u n ) , i = 0,...,n. (7.20) 

The resulting /3-functions indeed give rise to a NGFP with > 0, A* > 0, whose n-dependent 
position is shown in Table[2j Following the computation (|5.11j) - (|5.14p . we also obtain the stability 
coefficients of the NGFP. These are summarized in Table [3l In particular, we see that the NGFP 
attracts only the couplings uo, 1*1,1*2 as & ~~ * 00 > while it is UV repulsive for Ui,i = 3,..., 6. 
These results are in complete agreement with the analysis carried out in [22] and provides a nice 
confirmation of their findings! 16 ! 

8 Summary and discussion 

In this paper, we used the functional renormalization group equation (FRGE) for quantum grav- 
ity based on the effective average action to derive a renormalization group (RG) equation for 
/(-R)-gravity, where the gravitational Lagrangian is based on a arbitrary function of the scalar 
curvature. One prime feature of this flow equation is that its operator traces are invariant under 
constant rescalings of the scale-dependent functionals fk(R) — * cfk(R)- Based on this property, 
we were able to deduce that non-local interaction terms built from the curvature scalar can be 

16 The small numerical deviations between the results given in Table [5] and [3] and [25] can be traced back to the 
contribution of the single mode term in (|3.21[) . 
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consistently decoupled from the gravitational RG flow. In particular, the coupling constants ap- 
pearing in the cosmologically motivated ln(i?)- and i?~"-truncations (jl.4p and (jl.5p are constants 
along a RG trajectory, and thus the corresponding interactions are not created dynamically as a 
quantum gravity effect. 

Combined with the results [13j . these observations suggest that all gravitational interactions 
which do not occur in the heat-kernel expansion can be consistently decoupled from the gravi- 
tational RG flow. This does not imply, however, that is local at all scales k. In fact, more 
complicated non-local interactions built from the curvature scalar and the (inverse) Laplace op- 
erator D 2 , of the type / d d x^fgR(D 2 )- 1 R or, possibly, / d d Xy/gR\n{-D 2 )R, do appear in the 
late-time expansion of the heat-kernel [57] and are therefore most likely dynamically generated. 
Including such interactions in a truncation ansatz for 1^ is, however, beyond the scope of the 
formalism employed in this paper and will require its generalization to background metrics which 
are non-Einstein. Nevertheless, it would be very interesting to extend the present formalism 
in such a way to allow for the inclusion of such terms in the truncation ansatz, both from a 
phenomenological perspective (see, e.g., [59]), and from a more fundamental point of view, with 
respect to curing the conformal factor instability [491 150] or finding a dynamical solution to the 
cosmological constant problem. 

Motivated by recent attempts to explain the origin of dark matter from a modified theory of 
gravity which includes interactions that become dominant as the curvature scalar becomes small, 
the main focus of the present paper was the study of the non-perturbative RG flows arising 
from non-local extensions of the Einstein-Hilbert truncation containing J d d Xyfghx{R) (ln(i?)- 
truncation) or J d d x^J~gR~ n (i?~ n -truncation) terms in the ansatz for 1~V Including such non-local 
interactions in the truncation ansatz significantly changes the structure of the non-perturbative 
/3-functions obtained in the Einstein-Hilbert truncation. Particularly, the IR singularities of 
the RG trajectories leading to a positive cosmological constant (found in the Einstein-Hilbert 
truncation and including the "RG-trajectory realized by nature" [56]) are resolved. Notably, the 
i? -n -truncations with n > re cr it contain RG trajectories which are dragged into the NGFP for 
k — > oo, while giving rise to a positive cosmological constant in the IR as k — ► 0. 

Probably the most surprising result of this paper is the infrared fixed point (IRFP) appearing 
in the ln(i?)-truncation. For RG trajectories sitting on the plane of vanishing non-local coupling, 
this infrared attractor dynamically drives both Newton's constant and a positive cosmological 
constant to zero as k — > 0, abolishing the need of fine-tuning these couplings at the initial 
scale. While it is premature to claim that this provides a fullfledged solution to the cosmological 
constant problem, it is a first explicit example illustrating how could be dynamically driven 
to zero through strong quantum gravity RG effects in the infrared. In light of the particular 
characteristics associated with the inclusion of non-local interaction terms built from the curvature 
scalar only, it is crucial, however, to show the existence of such an IRFP at the level of more 
sophisticated truncations in order to establish the physical viability of this attractor. 

At this stage, we feel obliged to comment on the relation between the results obtained by 
studying the RG flows in this paper and the asymptotic safety conjecture for quantum gravity. 



46 



Notably, all our truncations give rise to a NGFP which is UV attractive for both Newton's 
constant and the cosmological constant. However, for some truncations this fixed point is located 
at negative values of Newton's constant, g* < 0. At first glance, this seems to be problematic 
for the asymptotic safety scenario. There, the NGFP is expected to arise at a positive value of 
Newton's constant, so that it provides a well-defined UV limit for the physically relevant RG 
trajectories with g > 0. In order to interpret our findings, it is then illustrative to compare the 
anomalous dimension r]M of Newton's constant obtained in the ln(iZ) and R~ n truncation with the 
one obtained in the Einstein-Hilbert truncation [4|. While 77 at in the Einstein-Hilbert case receives 
contributions from arbitrary powers of g, the t/tv arising in our non-local truncations are ony linear 
in g. Thus, it is likely that the non-local truncations considered in this paper are insufficient to 
correctly capture the RG flow of gravity in the UV, even though they are technically extensions 
of the Einstein-Hilbert truncation. In light of this, the shift of the position of the NGFP could 
very well reflect a "truncation artefact" , in the sense that having g* > requires a more refined 
truncation ansatz or the consistent decoupling of the new non-local interactions. Indeed, following 
the later procedure, the Taylor expansion of our flow equation for /(-R)-gravity reproduces the 
results of [22J, verifying the existence of the NGFP in truncations including interactions up to 
order R 6 . Ultimately, however, one would like to establish the existence of the NGFP (or a 
possible IRFP) from the partial differential equation governing the RG flow fk(R) by avoiding 
the expansion of fk(R) in a Taylor or Laurent series altogether. In any case, we feel that the 
exploration of the theory space underlying the FRGE approach to quantum gravity is still in its 
infancy and one might wonder which features still await discovery. 
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A Heat-kernel expansion on the ^-spheres 

In order to evaluate the traces on the RHS of the flow equation (|3.2ip , we make use of standard 
heat-kernel techniques. This appendix collects the relevant formulas. We thereby follow the 
discussion and conventions [7] extending their results by adapting some formulas from [60] . 

A.l Heat-kernel coefficients for unconstrained fields 

As the key formula for evaluating operator traces built from the covariant Laplacian —D 2 , we 
use the early-time expansion of the heat-kernel 



Tr 



d/2 

-it{D 2 +Q) ' ' " 



i \ ' f _ 

— j d d x^/g {ti a Q - ittr a 2 (x; Q) -t 2 tra 4 (x;Q) + ...} . (A.l) 
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Here, Q denotes a matrix potential, which we take to be proportional to the unit matrix in 
field space, and the a 2 k are the heat-kernel or Seeley coefficients!^! Furthermore, "tr" denotes a 
matrix trace in field space running over the tensor indices of the fields on which — D 2 acts. The 
derivation of the partial differential equation governing the RG flow of /(-R)-gravity in Section [7J 
then requires the knowledge of (some of) the a 2 k up to k = 4. 

The first two coefficients in the expansion (jA.ip are universal in the sense that they do not 
depend on whether — D 2 acts on unconstrained scalar, vector, or symmetric tensor fields: 



a 



a 2 (x;Q) = P , 



P--Q + -R1. 
b 



(A.2) 



Here, 1 denotes the unit matrix in the corresponding field space, i.e., 1( ) := 1, [l^)] 

fiu '■ — 9 put and 

[l(2)]pupa := 9pp9ua- Starting from a,4(x;Q), the heat-kernel coefficients include the commutator 
of two covariant derivatives on field space, 
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The corresponding terms then give different contributions to the a 2 k, depending on whether —D 2 
acts on scalars, vectors, or tensor fields. For a±(Q) we have (viz. [7]) 
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(A.4) 



for spin 0,1, and 2, respectively. In order to find the contributions from ae(x;Q), we start from 
the general formula [60] and restrict to spherically symmetric backgrounds. Taking into account 
the different normalizations this yields 
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For = 4 we need the scalar coefficient as(Q) only. Using the general formula [60J this becomes 
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Consulting the expansion of the traces (|7.1|) , (|7.2|) , and (|7.7|) , we then see that the knowledge of 
these heat-kernel coefficients suffices to evaluate the RHS of (|3.2ip in Section [71 

Note that the Heat-Kernel coefficients defined in [60] differ from the ones in (|XTj) by a factor (fe/2)!. 
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Table 4: Eigenvalues and degeneracies of —D 2 on the d-sphere. Here T lm (x), T l ™{x) and 
form a complete orthonormal basis for the — D 2 -eigenfunctions for scalars, transverse vectors and 
transverse traceless symmetric tensors. We refer to [7j for more details. 



A. 2 Heat-kernel coefficients for fields with differential constraints 

Now, the fields fiE v and are subject to the constraints given by the TT decomposition and 
thus, to compute their heat kernel coefficients, we must modify the formulas for spin 1 and 2 
fields accordingly. To do so, we first note that the TT decomposition f|2.Tj) and f|2.9j) implies the 
following relations between the heat-kernel traces over the constrained and unconstrained fields 
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Here we set Q = gJRl. There are two types of terms appearing on the RHS of these equations, 
complete traces with respect to the constraint fields and contributions from a discrete set of 
eigenvalues, which are required to complete the traces. The relations (|A.7|) then allow one to 
compute the heat-kernel coefficients for the constraint fields from their unconstrained counterparts 
given in the last subsection. Defining a2k\s '■= tr[a2ifc|J with s = 0, 1T,2T denoting the spin of 
the constraint field and performing the trace over vector indices explicitly we find 
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and 



«o| 2T = \( d - 2)(d + I) i 

(d-5)(d + l)(d + 2) 
a2 ' 2T " 12(^1) 

_ (d + l)(5d 4 - 22d 3 - 83d 2 - 392d - 228) 2 ^ A - 9 ^ 
a4,2T " 720d(d-l)2 R ' 

_ (d + l)(d + 2)(35d 5 - 287d 4 - 93d 3 - 7765d 2 + 1966d - 14016) 3 
° 6l2T " 90720 d 2 (d-1) 3 ^ ' 

for the heat-kernel coefficients of Tri^[e^ ltD2 ] and Tr2T[e _lt£,2 ]j respectively. Since the traces 
(|A.7j) are evaluated under a Fourier integral (cfg. eq. (I3.25|) ) it is also straightforward to deal 
with the discrete mode terms. Carrying out their (inverse) Fourier transform, these become 
proportional to <5-functions which allow for an easy evaluation of the subsequent s-integral. 

With regard to the derivation in Section [7] it is also convenient to give the numerical value of 
the heat-kernel coefficients for the special case d = 4. In the scalar sector, setting Q = 0, we have 

1 29 37 149 

a | =l, « 2 | =g^ «4| = ^^ 2 , «6lo = ^^ 3 , 6531840 ( A - 10 ) 

while ([Oil and ||OJ) yield 



1 67 

tr[ao| 1T ] =3, tr [a 2 | 1T ] = -R, tr [o 4 | 1T ] = r2 i 

5 271 7249 

tr[a | 2X ]=5, tr [a 2 | 2T ] = -- 12, tr [a 4 | 2T ] = -— 12 2 , tr [a 6 | 2T ] = r3 



(A.ll) 



Furthermore, when writing down the flow equations in sections and El we used the following 
shorthand notation for the d-dependent factors of the heat-kernel coefficients, obtained by setting 
12 = 1 in the formulas above 

C| fc (d) := a 2fc | (l? = 1) , CH{d) := a 2k \ 1T (R = 1) , C 2 fe T (d) := a 2fc | 2T (12 = 1) . (A.12) 



B Threshold functions 

This appendix collects the details about the standardized dimensionless threshold functions 
&n(w), &n(w) and T^ ;m (u, v, w), T p n Ju, v, w), which were used to encapsulate the cutoff-scheme 
dependence of our flow equations in the main part of this paper. We start by giving their general 
definitions and some of their properties in Appendix IB.H before specializing to the optimized 
cutoff [61 \ in Appendix IB. 2 1 



B.l General definitions 

Naturally, the /^-functions derived from the FRGE depend on the choice (of the scalar part) of 
the IR cutoff operators (|3"T2j) . R k = k 2 R^\- D 2 / k 2 ) . At the level of the flow equations, it 
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is convenient to encode this cutoff-scheme dependence in the standard dimensionless threshold 
functions 

(B.l) 



r(n)7 (y + J R( )(y) + «;) ? 

Here, i? (0) (y) 

is a dimensionless profile function determining the shape of the IR cutoff, while the 
prime denotes the derivative with respect to the argument. We note that both ^(w) and <&n(iu) 
satisfy a recursion relation when taking derivatives with respect to their argument, 

d w ^ n {w) = -pK +1 H , d w $P(w) = -pK + \w) . (B.2) 

In order to evaluate the last trace in (|3.2ip , a slight generalization of these threshold functions 
is required 



r ( n ) Jo (u (y + R(°) (y)) 2 + v(y + R<°) (y)) + w) ? 



(B.3) 



r ( n ) Jo (u (y + R(°Hy)) 2 + v(y + R^{y)) + wY 



These are similar to the threshold functions in [9], but adapted to the form of the traces appearing 
in (|3.2ip . The standard threshold functions (|B,1|) are contained in (|B.3h as a special case, 

T%(0, 1, w) = *p («,) , f p 0;0 (0, 1, «,) = & n (w) . (B.4) 

Obviously, the generalized threshold functions are homogeneous of degree — p under a common 
rescaling of their arguments 

TP m (A«, A«, Aw) = \- p T* m (u, «, «;) , T* ^(Au, A«, Aw) = A~ p T^ m >, u, w) . (B.5) 

Furthermore, they satisfy recursion relations similar to (|B.2|) 

^ T n,m( u ^> w ) = -pT^Ji f2 («,u,«j), 9 u T^ mZ (n,v,w) = -pT^^n.w), 
^ T n,m(^^^)= -PTSJ +I (u,»,!i;) , 9j[ mii (M,t),«)) = -pT^ +M (ti,!),w), (B.6) 
Sui T^ m («, u, w) = -pT^(u,D,w), ^" ^ jm)i (u,D,u)) = -pf^j^U,^) . 

The recursion relations (|B.2|) and (|B.6P are particularly useful when performing the series expan- 
sion of -R-dependent argument as, e.g., in Subsections 13.4} ISTT1 and [531 

B.2 Evaluation for the optimized cutoff 

For numerical evaluations, one needs to fix the form of the cutoff-scheme explicitly. Throughout 
this paper we thereby employed the so-called optimized cutoff [61] 

R k (z) = (k 2 - z)6{k 2 - z), z = -D 2 , (B.7) 
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with profile function 



R i0 \y) = (l-V)e(l-V), y = -D 2 /k 2 . (B.8) 

On the conceptual side, the main virtue of this cutoff is the minimization of the (modified) 
Ward-identities, thereby leading to more reliable and accurate RG flows. On the technical side, 
the optimized cutoff allows to carry out the integrals appearing in the (generalized) threshold 
functions analytically. Substituting (|B.8p into (|B.lj) . one obtains 

^ ) = T^(TT^' ^ w) = v^)jihy- (R9) 

Similarly, the generalized threshold functions (|B.3p become 

1 1 



Tin + 1) (u + v + w 



~ p (-1)" (5(-l,n,l + l) + (5{-l,n + l,l + l 



(B.10) 



[u,v,w) 



r(n) (u + v + w) 



Here, j3{— l,n, I) denotes the incomplete beta function. For fixed values n,l, these become con- 
stants. In particular we have /3(— 1,2,2) = 5/6 and /3(— 1,3,2) = —7/12. Observe that the 
optimized cutoff induces a degeneracy in the generalized threshold functions, in the sense that 
their RHS is m-independ. Making use of the properties of the incomplete beta-function, one can 
also verify explicitly that (|BT9|) and ([BlOl) agree m the limit (|B~4j) . 
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